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Abstract 

Roads to turbulence in open- flow shear layers are interpreted as 
sequences of often competing instabilities. These correspond to 
primary and higher-order restructurings of vorticity distributions 
which culminate in convected spatial disorder (with some spatial 
coherence on the scale of the shear layer) traditionally called 
turbulence. Attempts are made to interpret these phenomena in terms 
of concepts of convective and global instabilities on one hand, and of 
chaos and strange attractors on the other. The first is fruitful, and 
together with a review of mechanisms of receptivity provides a unifying 
approach to understanding and estimating transition to turbulence. In 
contrast, current evidence indicates that concepts of chaos are 
unlikely to help in predicting transition in open-flow systems. 
Furthermore, a distinction should apparently be made between temporal 
chaos and the convected spatial disorder of turbulence past Reynolds 
numbers where boundary layers and separated shear layers are formed. 
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1. INTRODUCTION 

Information on instabilities and transition to turbulence in 
shear layers has amplified exponentially in recent years and useful 
earlier summaries Morkovin (1978), Reshotko (1976, 1984), Arnal (1984), 
Ho and Huerre (1984), Schlichting (1979) and others can serve only as 
essential points of departure for a novice in the field. 

This report represents a selective account of the more recent 
insights into transition mechanisms in open- flow systems . Fig. 1 . as 
well as some new questions and disagreements that have arisen and that 
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may challenge innovative researchers. This up-to-date perspective is 
presented by focusing on vorticity fields: the transition process 
consists of a sequence of often competing instabilities which 
correspond to primary and higher-order restructuring of vorticity 
distributions. Use of composite figures was chosen as a key tool to 
convey and summarize overall features and contrasting developments. A 
careful student of the subject will supplement these views with others, 
especially more mathematical approaches such as that of Bayly, Orszag 
and Herbert (1988) . Spatial constraints limited the number of 
references; mostly they were chosen as particularly illustrative and 
suggestive. From these references and discussions in them a broader 
bibliography can be readily compiled. 

Two recurrent themes wind their way through the text: (1) 
response of shear layers to external stimulation, "natural” (indigenous 
to their specific environments) and intentional, by vibrations, sound, 
protuberances, grids, etc. for controlled modifications of behavior, 
and (2) comparison of instability and turbulence behavior in open- flow 
and closed- flow systems. The first theme rests on the recent 
clarification of concepts of convective and global instability and 
proves itself very fruitful. It is coupled with the mechanisms of 
receptivity by which the external disturbances get internalized as 
unstable intrinsic shear-layer vorticity modes. The second theme aims 
at the possible applicability of the concepts of chaos and strange 
attractors to the open- flow systems. A simple introduction to these 
concepts is given in the Appendix with emphasis on the first theme: 
response to harmonic forcing. 

Broadly speaking, open- flow systems, Fig. 1 . have two 
distinguishing behavior patterns not seen in the inherent oscillatory 
motions of non- fluid systems: downstream propagating primary 
instabilities within a broad band of frequencies , and (at the end of a 
sequence of higher-order instabilities) spatially disorderly turbulence 
with some spatial coherence on the scale of the shear layer (the 
coherent "large scale" structures) . All convectively unstable open 
flows at higher Reynolds numbers display their particular 
disturbance -conditioned sequence of instabilities successively in x 
(the streamwise coordinate) at the same time t. Restructuring of 
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vorticity probably continues in all shear layers with x-varying 
thickness. Only far downstream in internal flows, Fig, lb , can a true 
equilibrium of turbulence possibly take place. 

In contrast, within closed- flow systems the flow interacts with 
itself at all times . Under this iterative self - interaction, 
equilibrium states at constant Reynolds number (including turbulent 
states) are generally reached in a short time at all spatial 
locations. A fixed Reynolds or Rayleigh number and a fixed stable 
oscillatory equilibrium can usually be set at will and maintained in 
time by the experimenter-designer. It will be seen that these 
differences impair the applicability of the concepts of chaos to 
open- flow systems. 

The interweaving of the two themes and the structure of the paper 
can be gleaned from the Table of Contents and the sequence of figures. 
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2. ROADS TO TURBULENCE IN OPEN-FLOW SYSTEMS 

(2.1) Dehomogenizing , linearized primary instabilities . 

Instabilities arise primarily in thin shear layers 
generated by viscosity in flows past solid walls: boundary 
layers, wakes, jets, mixing layers, pipe and duct flows, etc. 

At lower Reynolds numbers, Re, the laminar velocity profiles 
(such as those in Fig . 2 for mixing layers, flat and concave 
boundary layers) form quasi- two-dimensional continuous vorticity 
distributions which are nearly homogeneous in the streamwise x 
and spanwise z directions. As Re increases, vorticity dynamics 
lead to locally increased concentrations Aw associated with the 
first or primary instability as sketched in mid Fig, 2 , This 
first restructuring of vorticity consists thus of dehomogenization or 
lumpization of the w field in which total vorticity is preserved. 

In the case of the concave boundary layer, incremental positive 
and negative streamwise vortices, named after Goertler and Taylor, 
grow in x, Fi g , 2 right, while the flow remains steady. In the case 
of the flat-plate boundary layer, positive and negative intensifi- 
cations of spanwise vorticity Aw are convected downstream at speed c. 
The dehomogenization together with convection render the flow unsteady. 
However, an observer traveling downstream with speed c would see a 
steady growth of the concentration in the so-called Kelvin cat's eyes, 
located slightly above the critical height y^(for which U(y cj J- c) , 
center of Fig . 2 . In the case of the antisymmetric mixing layer on the 
left of Fig . 2 spanwise Aw grow in time and remain in place. 

(2.2) Amplification and damping of first and higher normal modes 

This restructuring in open- flow systems generally proceeds from 
vorticity disturbances (nonhomogeneous in space and/or time) induced by 
small perturbations in the environment of the given shear layer. The 
initial growth rate of the vorticity concentrations can be computed 
from the linearized Navier- Stokes equations for the local geometry of 
the layer. When excitation is purely sinusoidal . the response consists 
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of a superposition of normal or eigen modes associated with the local, 
nearly parallel structure of the shear layer. In general, all but the 
lowest mode are damped. Above a geometry- dependent minimum critical 

Reynolds number, Re ^ , the lowest mode is selectively amplified in 
relatively narrow frequency and/or wave number band. The corresponding 
y-distribution of vorticity, velocity components, pressure, etc are 
called eigen or characteristic functions and can at times be verified 
experimentally. They grow proportionately to the cosine- exponential 
combinations shown under the sketches in Fig. 2 . Thus, the incremental 
vortices in the mixing layer case intensify in place at the temporal 
exponential rate q^c^ per unit time (where wave number is the 
inverse of the wavelength A multiplied by 211 and c^, c^ are the real 
and imaginery components of the complex "propagation speed "c). For 
the other two cases in Fig. 2 . the exponential amplification rate is 
spatial in character, -cr^ per unit length, where a r » a £ are 
the real and imaginary parts of the complex wave number a, again 
computed as eigenvalues for the given system of linear equations. 

(2.3) Nonlinear developments: secondary and higher instabilities 

It has been recently established that nonlinear effects, neglected 
in the above computations, can initiate secondary resonance mechanisms 
when the r.m.s. x - velocity fluctuations u' reaches relatively low 
levels of 0.3 percent of the mean velocity at the edge of a 
flat-plate boundary layer. In mixing layers secondary restructuring 
commonly begins at primary u' levels near 0.001* U , while the primary 
vortex growth and roll-up continues nearly exponentially to saturation 
levels on the order of 0.1*Ue. (The cited u' levels convey relative 
not absolute energetics. In experiments, even numerical ones, the 
energies are contained in nonstandarized finite band-widths, which 
vary among experimenters.) It is important to note that slower, hardly 
diagnosable three-dimensional tertiary restructuring may commence 
almost hand in hand, with the secondary vortex-pairing instability. 

These observed cases thus do not conform to the scenario postulated 
by most mathematicians, namely that the primary instability should 
first reach a nonlinear equilibrium . Such a "saturated" equilibrium 
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flow would form a new base flow which, in principle, could be perturbed 
and linearized in order to characterize the secondary instability modes 
bifurcating from the equilibrium flow. One reason for the growth of 
secondary instabilities from non- equilibrium states is the fact that 
these shear layers constantly grow in thickness. That thickness is 
taken as constant in the bifurcation theory. Another reason may be 
that in the strictly parallel boundary layer case the bifurcation is 
"subcritical" . Subcriticality implies that the finite -amplitude 
equilibrium, if achievable, would be automatically unstable. Herbert 
(1988) describes in great detail the experimental and theoretical 
research which resulted in the current understanding of secondary 
instabilities in boundary layers (summarized in the center of our 
Fig 15 . ) Even though the instabilities start from non-equilibrated new 
base flows, the assumption of a saturated primary leads to a very 
satisfactory agreement with experiments, Herbert (1988). 

(2.4 ) Parametric instabilities 

Herbert also stresses that all the first restructurings bring about 
systems of equations in which some coefficients of the unknowns vary 
periodically in space or time, say with period P, as sketched in Fig. 2 . 
The resulting linear Floquet equations exhibits powerful parametric 
instabilites , with periods (m/n) P, where m and n are integers. Thus 
there are many possible competing secondary instabilities, depending on 
initial conditions and on the "saturated" amplitude of the primary mode 
(threshold effects). Consequently, for each physical realization, such 
as those in Fig. 2 . this unified mathematics will be reflected in 
different competing geometrical vorticity formations . In mixing 
layers, the especially strong subharmonic instability with m - 2 and n 
n - 1, may correspond to pairing of two-dimensional primary vortices 
with total vorticity conserved or to a three-dimensional corrugation 
of the primary vortices, with the quintessential new-vorticity 
generating mechanism by "tilting and stretching"; see theory of 
Pierrehumbert and Widnall (1982) and experiments described 
by Ho and Huerre (1984). For the flat plate boundary layer in Fig: 2 . 
the three-dimensional secondary instability (m - 2, n - 1) is present 
in the Craik (small amplitude) and Herbert (finite amplitude) 
mechanisms. The earlier explored secondary instability, the K 
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mechanism, studied experimentally by Klebanoff, Tidstrom and Sargent 
(1962), corresponds to m - 1 and n - 1, and to a higher threshold 
amplitude of the primary instability. 

(2.5) When does a turbulent state set in? What is it? 

As indicated in the summary Fig. 3 . higher instabilities, i.e., 
rapid restructurings of vorticity are observed as Re increases, often 
simultaneously on different scales. In closed flow systems, such as 
flows between rotating cylinders, the successive instabilities are 
relatively clear cut. Yet there is competition between the many 
possible configurations. Also, localized regions of turbulence occur, 
especially in counter- rotation, see Andereck, Liu and Swinney (1986). 
In open- flow systems with shear layer thickening in x or t, lack of 
resolving power of our instruments limits identification of 
instabilities beyond the tertiary, - if any are present. In a short 
time or distance past the occurrence of the tertiary we are generally 
unable to distinguish the state of the flow from a classically 
turbulent one. According to Tennekes and Lumley (1987) and R. W. 
Stewart (1968) four key syndromes characterize turbulence encountered 
in engineering, Fig. 3 . : (1) irregularity, (2) three-dimensional 
vorticity, (3) diffusive motions far in excess of molecular mixing, and 
(4) dissipation. The first three are useful in diagnostics: they imply 
broad band (three dimensional) spectra. In motion-picture 
visualizations the onset of turbulence at higher Reynolds numbers is 
commonly seen as a sporadic sudden localized increase in diffusion 
which oblitarates small regions of laminar streak lines as an 
" incipient turbulent 3D spot " is formed in the final stage of 
instability . 

Such emerging turbulent 3D spots need to be distinguished clearly 
from the turbulent spots and patches studied experimentally in boundary 
layers and ducts which are often called Emmons' spots. Emmons' spots 
grow in an environment of nearly quiescent laminar layer whereas the 
neighborhood of the incipient turbulent 3D spots, though still laminar, 
is highly unsteady and contains significant three-dimensional 
vorticity. Dynamics of their subsequent growth are unlikely to be 
comparable . 
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To what extent are the classically turbulent states in open-flow 
systems (as characterized by the above four syndromes) equivalent to 
" chaotic " states characterized by " strange attractors " in the modern 
theory of dynamical systems? What deeper insights can the theories of 
chaos bring to our dealing with turbulence in open flow systems and in 
particular to our engineering task of estimating the x- location of 
transition to turbulence? The 1988 Spring issue of a certain alumni 
magazine quoted an eager researcher at the university: "...we are now 
able to predict the point at which a laminar flow will start to become 
turbulent..." and "This applies to the turbulent flow encountered by 
any type of vehicle in any fluid". We note that the four- syndrome 
characterization of turbulence refers to localized regions of the shear 
layer and allows for turbulence to coexist in proximity with laminar 
flow as in both types of turbulent spots above. In the theory of 
dynamical systems the designation chaotic refers to the state of the 
whole system - the whole shear layer. Locating the x- onset "point" of 
turbulence in a vehicle boundary layer requires additional assumptions 
which should be rigorous. Swinney and Gollub (1986), who pioneered the 
experimental study of chaos in closed- flow systems, aware of such 
difficulties, state in a review article: "It is not yet clear whether 
the strange-attractor concept is useful for open systems..." 

Before attempting to address the questions in the preceding 
paragraph, we need to clarify in Section 3 what is meant by convective, 
absolute and global instability in open- flow systems. A simple 
illustration of the concept of strange attractor and some philosophical 
implications of chaotic phenomena is found in the Appendix. 

(2.6) Re Re .Re, .and bvt>ass transition 

Two more physical concepts which distinguish between intrinsically 
laminar flows and decaying turbulent flows are signaled in Fig . 3 . 

There is ample evidence, Morkovin (1984), that boundary layers below a 
Reynolds number, R e m ^ llL * cannot sustain intrinsic self -energizing 
turbulent motions on the scale of the local shear layer thickness 
6 (x) . In turbulent boundary layers, the unsteady energy dissipated 
near the wall is restored through intermittent "bursts" which produce 
negative Reynolds stresses - u\7 and recoup energy from the mean 
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flow. Two types of events associated with bursts, ejections and sweeps, 
are identified in Fig. 3 . The ejections in particular have the 
character of local instability which evidently is absent below Re m £ nt 

Any locally forced turbulence, say by a spark or local separation, 
decays rapidly. Even when vigorous free-stream turbulence is convected 
toward the boundary layer as in Fig la . the energizing "bursts” (with 
scales smaller than 6 ) are not excited. Rather, the boundary layer 
experiences some three-dimensional motions on scales larger than <5 and 
some mild increases in mean wall stress and heat transfer rates . 
However, the boundary layer remains intrinsically laminar. This 
randomly unsteady quality is conveyed by the term " buffeted laminar 
layer " in Figs, la and 3. For the given shear- layer geometry the 
earliest possible transition, important in design, occurs at Re m £ nt 

The observed minimum Re for self - sustained turbulence in four types 

of flows are compared to corresponding critical Re for first laminar 

amplification of inf initesinal disturbances in Fig. 4 . In presence 

of moderate to small environmental disturbances of turbulence, Re_ in 

tr 

boundary layers is generally observed downstream of Re cr > at distances 

corresponding to inferred total amplifications of roughly 100 to 10,000 

of the linearly most dangerous disturbance. The observed approximate 

coincidence of Re . and Re of about 200 (based on momentum 
mint cr 

thickness 6 ) for the flat plate has remained unexplained for forty 
years. However, the expected sequence Re m ^ nt < Re cr < ^ e tr S enera Hy 
observed in boundary layers on bodies with finite thickness when 
disturbances are moderate. 

Occasionally, Ra tr does not appear to be related to Ra cr and may 
even occur upstream of Re cr - Obviously, such cases can play havoc with 
design, as it did in the humbling case of the heat-sink nose on the 
first reentry vehicles in 1956, Murphy and Rubesin (1966), Morkovin 
(1984). (This is called the blunt-body paradox in Fig. 15 .) On blunt 
bodies at high speeds the transition path to turbulence bypasses to 
this day all known instability theories and documented mechanisms 
(although highly accelerated thin boundary layers are known to be very 
sensitive to surface roughness). We can expect large disturbances 
(local wall distortions, vigorous oncoming turbulence) to cause bvpass 
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transition as outlined in Fig 11 . However, the mildly disturbed 
"naturally occurring" bypass transitions Re m ^ nt < ^ e tr < Re cr ’ 
pipes and ducts in Fig. 4 also remain unexplained. These flows are more 
subcritical to axisymmetric and two-dimensional disturbances than 
boundary layers (Section 2.3), but experiments point to 
three-dimensional disturbances and more subtle causation. 

(2.7) Relaminarization 

The final important concept in the summary Fig. 3 deals with 
relaminarization, which has been extensively documented and analysed by 
Narasimha and Sreenivasan (1979). One of the several relaminarization 
mechanisms they describe operates in highly accelerated turbulent 
boundary layers. The favorable pressure gradient increases dissipation 
near the wall and stabilizes the inner layer. The turbulence -re- 
energizing bursts diminish in frequency and stop altogether. A thin 
new laminar layer grows from the wall within the old boundary layer, 
which now decays and becomes wake-like. The outer-layer turbulence 
cannot stop instantaneously; it is merely "starved". In fact, the 
growing inner layer constitutes a laminar layer "buffeted" by the 
decaying turbulence of the outer layer. It remains so until it grows 
past its Re mint » Sternberg (1954), or until the favorable pressure 
gradient is relaxed and bursts recommence. Again, the wall stress and 
the heat transfer rate in such relaminarized buffeted inner layers 
appear to be less than those expected in turbulent layers, justifying 
the term "locally relaminarized". 
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3. CONVECTIVE INSTABILITY AND STREAMWISE RESONANCE 
(3.1) Local convective and absolute instability 

The normal -mode responses to the specialized harmonic excitation of 
the linearized primary instability, Fig. 2 . were termed spatial or 
temporal . A more fundamental instability response of a shear layer is 
that to an impulsive forcing disturbance at t ■ o, localized at a the 
origin of the streamwise and spanwise coordinates. In terms of Dirac's 
impulse functions, response is desired to the forcing term 6 ( t) 6 (x) 6 (z) 
on the right side of the otherwise homogeneous linearized stability 
equation. This search for a Green's function of the linearized system 
represents an idealization of the experiment of Gaster and Grant 
(1975). By repetitive excitation and judicious sampling these 
researchers measured the unsteady response of a Blasius boundary layer 
to short weak pressure puffs from a very small hole in a flat plate. 

At any time after the puff, the boundary layer was disturbed only over 
a small x-z area which grew slowly in size as the disturbance field 
propagated downstream. In accordance with the filter-amplifier 
characteristics of the layer, the disturbance took the form of x-z wave 
packets which grew at first nearly exponentially in amplitude within 
the disturbed, propagating “ laminar spot " . At later times nonlinearity 
and higher instabilities modify the nearly exponential linearized 
primary instability behavior. Such localized exponential growths 
within downstream propagating spots are defined as being convectivelv 
unstable . 


Figure 5a . based on Huerre and Monkewitz (1985), and Huerre (1987) 
illustrates the linearized convective instability in one less dimension 
than the Gaster-Grant experiment. As in Fig. 2 the response is 
exponential, but in both x and t simultaneously, yielding a growing 
wave packet within the propagation wedge shown in Fig. 5a . Outside the 
wedge the factor exp(u>^t) with negative damps out the response. The 
wave pocket transmits disturbance energy at the group velocity c 


du> 


r /5k^ which is associated with a ray x-c^t within the wedge 


g 
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In some flows, illustrated in Fig 5b . the group velocity becomes 
zero and the disturbance energy remains centered on the point of 
impulse, x - o. The disturbances then grow exponentially in place and 
disperse somewhat in x according to the normal-mode dispersion 
relation. Such linearized flows are called absolutely unstable . In 
open-flow systems, most of the mean vorticity of the shear layer 
continues to be convected downstream but the energy, associated with 
the vorticity disturbance field may be trapped, "working upon itself. 
We note paranthetically that this self - interactive character of the 
disturbance field is partially similar to the interactive conditions of 
Section A-4 for Couette and Benard-Rayleigh flows, which ultimately 
lead to strange -attractor behavior. Those two flows are indeed 
absolutely unstable, but in addition they are closed with respect to 
the mean vorticity, an important feature when -turbulence and phase 
jitter set in. 

(3.2) Linear resonance in x. nonlinear effects and pressure feedback 

The spot- like character of the impulse response requires both the 
circular frequency w and the wave number k (often denoted by a) to be 
complex in the spatio - temporal factor expi(kx-wt) in Fig . 5 . A 
functional relation w(k) (the dispersive relation in terms of the flow 
parameters) is imposed by linearized equations for the disturbance in 
the given flow. The absolute instability condition dw/dk - 0 of zero 
group velocity generally corresponds to a double root for the 
eigenvalues of w. Two disturbance vorticity waves propagate then at 
the same speed w^/ k r and potentially could "resonate". Koch (1985) 
postulated that vanishing of the u> derivative in the wake of a bluff 
body indeed corresponds to a resonance leading to the extra strength 
and regularity of the observed Karman-Benard vortex street. The 
resonant-like nonlinear regularity of the Karman-Benard street is now 
believed to be associated with global instability of the flow, 

Monkewitz (1988), for which the initial absolute instability is only a 
necessary but not a sufficient condition , see Section 3.4. 

The important definitions of convective and absolute instabilities 
are strictly valid only for small -disturbance linearized flows which 


-14- 


(3.3) 


are independent of x. The extra qualifying adjective " local ” for these 
instabilities underscores the suppression of the x dependence of the 
layer in the analysis, which then permits separation of variables and 
reduces the problem to the search for eigenvalues in y. Wakes of bluff 
bodies and jets with variable density (the other open- flow shear layer 
with absolute instability, Monkewitz and Sohn, 1986) unfortunately 
exhibit rather strong x-dependence and depart significantly from 
infinite parallel layers. Furthermore, the suspected resonance involve 
wave interaction in x with some unspecified reflection at the body and 
possibly elsewhere. 

In the bluff -body cases, the local instability criterion disregards 
completely the strong destabilizing adverse pressure gradient near the 
base which is suspected to contribute significantly to the magnitude of 
the pressure fluctuations. It is known that the common dominant 
low-frequency pressure disturbances in wind tunnels stem from irregular 
separation in diffusors with adverse pressure gradients. Instantaneous 
separation lines often tend to be nonuniform spanwise due to presence 
of u? x at moderate Reynolds numbers. However, resonances in x in shear 
flows tend to straighten out separation lines and lead to still larger 
fluctuations. There is nonlinear hydroacoustic feedback to upstream 
regions which reinforces or modifies any interactions (and feedback 
tuning) between vorticity waves. 

In other words, the shear layer we observe may differ substantially 
from the hypothetical ones in which infinitesimal disturbances 
presumably initiated the resonant processes. For instance, at Re^l20 
the streamwise velocity fluctation u' in a cylinder wake exhibits the 
nonlinear level of 0.0211^ right at the separation of the shear layer 
from the body, Nishioka and Sato (1978, Fig. 5). The pressure 
fluctuations seeded by the vigorous nonlinear vortex formations (see 
equations in Fig . 7 ) actually reach upstream of the cylinder and 
precondition the boundary layer on the cylinder upstream of the 
oscillating separation point. This exceeds any upstream effects which 
a vorticity wave interaction- resonance could have caused at the initial 
linearized level and is therefore unquestionably hydroacoustic 
in character. The same Nishioka- Sato figure documents that this case 
of local absolute instability (characterized by the purely temporal 
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initial instability in Fig. 5b ) evolves into a global instability in the 
form of a rapid purely spatial u' amplification in x for x/D up to 
approximately 2.5 and a slow decay thereafter. The maximum fluctuation 
u' (in excess of 0.31^) is located just downstream of the closure of 
the mean-wake recirculation pocket. This provides support for at least 
a contributory role of vorticity feedback in the final state of 
cylinder-wake resonance at low Reynolds numbers. 

The resonant bluff wake has an additional property, useful 
in diagnostics. It is robust with respect environmental disturbances. 
In other words, the nonlinear feedback system filters out external 
perturbations to the point that the Re-dependent frequency of the 
vortex street has been used successfully for velocity measurements at 
low speeds where Pitot tubes are inaccurate, Roshko (1955). The 
resonant near wake at low Re seems to have just two effective degrees 
of freedom and behaves like a second-order Van der Pol oscillator with 
a robust limit cycle described in Section A-2 and Fig A-la.b . 

(3.3) Noise response of convectivelv unstable layers 

Figure 6 summarizes the progression of idealizations leading to a 
model of nonlinear global instability, presumably mimicking the 
resonance of the cylinder near wake. First, for the linearized 
convectively unstable shear layer - present in the majority of open 
flows - the local parallel approximation exacts no significant 
penalties. Second-order nonparallel theories (when carried out 
correctly) provide relatively small correctons for Re^ and 
amplification rates. The convective character of the instability 
remains unchanged. 

Thus by definition of convective instability each successive 
disturbance "received by the layer from upstream" generates its own 
laminar spot which propagates downstream as it grows exponentially. In 
wind tunnel environments with continuous randomly modulated 
disturbances laminar spots tend to superpose while maintaining some 
random phase and amplitude nonhomogeneity in x and z. Sampled at any 
frequency within the amplified spectrum, the response resembles that of 
an averaged quasi- two -dimensionally spatially growing unstable wave; 
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see lower half of Fig. 13 . As these disturbances approach nonlinear 
levels, the amplified, noise- sustained structures (in the language of 
Deissler and Kaneko, Section A-4) are ready for secondary and higher 
instabilities, fed in part by the amplitude inhomogeneities in x and z. 

The convectivelv unstable shear layers lack any robustness with 

respect to environmental disturbances - on the contrary they are driven 

by them, probably even quite some distance into the fully turbulent 

region. The randomness and nonhomogeneity of the response is 

reflected in an extended region ARe^ within which the local onset of 

turbulence may take place. This is illustrated in the lower half of 
Fig. 13, where Re^ ^ e g refers to the upstream position of the onset 

and en ^ to the location past which all activity near the wall is 

burst-like, Section 2.6. Measured ratios ARe_ to mean Re in wind 

tr tr 

tunnels are often on the order of 0.4, but tend to be smaller in 
adverse pressure gradients. 

The vision of transition prediction bv the researcher quoted in 
Section 2.5 is thus strongly contraindicated for convectivelv unstable 
boundary layers "on any vehicle". 

(3.4) Absolute versus global instability 

Local absolute instability rests on the character of the eigenvalue 
problem in y. The occurrence of zero group velocity in the x-dependent 
shear layer signals merely the possibility of a resonance in x. The 
visual evidence of Thorpe (1971, Figs, 1, 5 and 7) for the anti- 
symmetric mixing layer demonstrates a lack of robustness in his 
absolutely unstable flow. Unquestionably there is no resonance. The 
layer responds irregularly to the initial nonhomogeneities, which are 
amplified into "noise - sustained" saturated structures, and proceed to 
irregular pairing and sporadic onset of turbulence. It would almost 
seem that the very lack of x-dependence in the mean mixing layer (in 
contrast to that of the bluff near wake) removed "the opportunity" for 
resonance, such as partial or full reflections of the vorticity waves. 

The situation seems to have been clarified on a partial 
differential equation model (the Ginzburg -Landau equation in x and t, 
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with a convective term UdA/3x, a cubic nonlinearity and one coefficient 
slowly varying in x) by Chomaz, Huerre, and Redekopp (1988). The 
linearized equation admits an undamped eigenfunction in x - a robust 
global instability - only when the x-extent of the equivalent local 
absolute instability reaches a finite critical size . Numerical 
solutions of the full nonlinear equation for impulsive white-noise 
conditions exhibited the successive behavior for the stable, 
convectively unstable, absolutely unstable, and globally unstable 
regimes as the parameter of the system was varied. 

Monkewitz (1988) in his study of local viscous linearized stability 
of a family of wake profiles associated the onset of the three unstable 
stages with Re^ of approximately 5, 25, and 47, respectivey. He took 
the empirical result for the spontaneous onset of the oscillations at 
Re D^ ^ aS criterion f° r the onset of global instability, a 

criterion not obtainable in parallel approximations. The common 
absence of oscillations in the convective regime 5 < Re^ < 25 is 
ascribed to the absence of environmental disturbances in the narrow 
excited range of frequencies. Taneda's 1963 visualizations of wakes 
downstream of slightly vibrating cylinders supports this 
interpretation; see discussion of Fig. 5 in Morkovin (1964). 

Judging by Fig. 16 of Nishioka and Sato (1974) for the 
recirculating wake, Monkewitz local conditions for absolute instability 
at Re^^v 25 correspond to the closure of the steady recirculation 
pocket at x^l.7D. Just below the spontaneous street onset at Re ^of 
47 the pocket extends to x/%*3.1D and presumably provides the 
increased critical range needed for the resonance in x. Just as in 
the model of Chomaz et al (1988) in this regime of absolute instability 
there is little qualitative distinction of behavior from that in the 
convective instability regime: external forcing is required for the 
instability to appear and grow. For instance minute cross-flow 
oscillations of the cylinder engender an initially linearizable 
exponential amplification in x for Re^ of 30 and 40, as documented in x 
in Figs. 6 and 7 of Nishioka and Sato (1978). 

Altogether, the picture now appears rather consistent despite the 
ambiguities, outlined in Section 3.2, associated with the definition of 
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local infinitesimal instability, the usage of mean profiles in 
highly disturbed flows and the role of the nonlinear hydroacoustic 
feedback. Particularly convincing were papers which studied the time 
evolution of the global instability after it was suddenly "switched on" 
and which verified that the whole field grew temporally with the single 
factor exp(w^ o t) in the linear regime, as one would expect from Fig. 5b . 
Experimental analytical verification came from Sreenivasan, Strykowski 
and Olinger 1987 and Provensal , Mathis, and Boyer (1987), while a 
numerical study for a thick flat plate with a square base was described 
by Hannemann and Oertel (1987), based on Hannemann's dissertation. 

The numerical experiments permit detail examination of some of the 
aforementioned ambiguities. In accordance with the scenario of Chomaz 
et al (1988), Hannemann (1988, Fig. 42) finds the flow damped at Re^ - 
80, where the range of local absolute instability extends from the base 
to x - 2. 3D, a range evidently insufficient to build up an amplifying 
eigenfunction in x. For Re^- 200, the disturbances in a 
quasistationary symmetric Navier-Stokes field grow everywhere by a 
factor or 10^ in accordance with the linearized factor expw. t while 
the free stream covers a distance of about 240D. The subsequent 
nonlinear evolution, starting at u' 0.007*11^ is marked by (1) gradual 
shortening of the mean recirculation pocket by eddy entrainment (Fig. 
20), (2) change in mean velocity profiles, (3) a 10% increase in the 

Strouhal frequency (Fig. 21), (4) evidence of hydroacoustic feedback to 
the boundary layer upstream of the base, (Fig. 28), and (5) exp(2u iQ t) 
growth of a varicose instability component of twice the frequency of 
the initially pure sinuous instability and movement of the unsteady 
vorticity maxima off the center line, (Figs. 24 and 19). 

(3.5) Convective -hydroacoustic global instability : numerical 
simulations 


In closing this section, we should note that a whole class of 
global instabilities of free shear layers is sustained by hydroacoustic 
feedback when they impinge on a solid boundary; see bottom of Fig . 6 and 
the overview by Rockwell and Naudascher (1979). When unsteady 
vorticity of a convectively unstable free shear layer impinges on a 
wedge, it generates pressure feedback to the point of separation of the 
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layer. Receiptivitv to unsteady pressures gradients at separation 
being very high, Morkovin and Paranjape (1971), significant unsteady 
vorticity is generated at the origin of the layer. This newly-born 
unsteady vorticity is amplified by the unstable shear layer and 
convected toward the wedge, thus completing the resonance loop. A 
singing teakettle commonly reminds us shrilly that such loops can build 
up substantial energy levels. There is resonance and at least one 
undamped eigenvalue in x; there is robustness with respect disturbances 
(except for conditions of competition between two neighboring 
eigenvalues) . 

The nonlinear generation of unsteady pressure is absent when 
Hannemann's global bluff -wake instability commences. But his Figure 28 
indicates that the upstream effects of the pressure grow with 
nonlinearity and that they probably contribute to vorticity control at 
the very origin of the wake in the later stages of the global 
instability at low Reynolds numbers discussed here. At high Re values 
the threshold for pressure feedback maybe very low. 

It is worth noting that the upstream traveling vorticity or 
pressure wave can be decaying in the negative x direction and yet 
sustain the global resonance. In the hydroacoustic case the pressure 
fluctuations decay in all directions from their origin, the exact 
manner depending on the geometry and the ratio of the total propagation 
distance to pressure wavelength (near-field versus far-field behavior). 
However, high receptivity followed immediately by large amplification 
in the layer and rather efficient conversion into pressure fluctuations 
at impingement, Rockwell (1983), provide sufficient regeneration and 
phase coherence to maintain the loop against more incoherent 
disturbances. We note that in the near-wake and variable -density- j et 
cases, the nature and dynamics of the resonant vorticity waves have not 
yet been identified theoretically or experimentally. The dominance of 
the spatial development, evident in the saturated state of the 
instability (Hannemann, 1988, Fig. 28; Nishioka and Sato, 1978, Figs. 4 
and 5, Unal and Rockwell, 1988a, Fig. 5) suggests that any upstream 
traveling vorticity mode may be dwarfed in the interactions and may 
possibly be decaying towards the separation point. 
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Paradoxically, accuracy of numerical experiments with strearawise 
amplification of free shear layers is endangered by inadvertent 
pressure feedback from the unavoidable conditions at the outflow 
boundary, OB. Ideally, differential relations could be imposed at OB 
by dropping the “viscous" highest order derivatives from Navier-Stokes 
equations. This would make OB transparent to oncoming vortical 
structures except for a thin unsteady viscous "shock-boundary 
adjustment layer" just upstream of OB. (This BC scheme has convergence 
and expense problems). Otherwise, constraints imposed at OB convert 
the convected dynamic field more or less efficiently to pressure 
fluctuations as in the case of impinging shear layers. Pressure 
feedback to the receptive separation region of the layer and 
subsequent amplification- convection provide an additional disturbance 
loop imposed upon the "natural" dynamics of the layer under scrutiny. 
From the com- puted results alone, it is difficult if not impossible to 
judge how serious the likely modifications of the dynamics may be - 
changes in phase, amplitude, frequency, wavelength, onset of higher 
instabilites , etc . 

Simpler differential conditions at OB with some promise are 

{|e - D Cy>- 1^} v - ° or ||e * u, |j} 2 - 0 

where U^is a convection velocity of large eddies near OB and V the 

velocity vector. They correspond to Taylor's hypothesis in hot-wire 

anemometry and presumably decrease conversion of the energetic eddies 

2 2 

into pressure fluctuations. The condition d u/dx - 0 used by 
Hannemann is adequate for low- amplitude vorticity waves in attached 
boundary layers with low receptivity, but its applicability to free 
shear layers with separation should be rigorously investigated. A 
similar problem, though probably less serious arises at the inflow 
boundary, IB. The practice of prescribing specified y,t or y,z,t 
velocity variations on IB may simulate tolerably the vorticity or 
pressure waves which penetrate the field of computations from upstream. 
However, since these prescribed values do not allow for any 
modifications by pressure fields spreading from within the computing 
domain or from OB, the condition automatically imposes a complicated 
velocity mode at IB, tantamount to special spurious reflections. 


I 
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4. NATURE AND ROLE OF RECEPTIVITY 

.(Axil Recent ivity- conditioned wea k global instability in mixing, layers 

In Section 3.5, we saw how a convectively unstable free shear layer 
can be converted into a globally unstable system by inserting into it a 
solid object which generates unsteady pressure feedback to the 
separation region. In its initial laminar segment many a mixing layer 
has the character of a weak global unstable system without intervention 
of an intruding body: the vigorous vortex pairing driving its 

secondary instability replaces the impingement region as the pressure 
source. Figure 7 summarizes the many coupled processes involved in the 
complicated two - frequency feedback loop as documented by investigators 
at the Illinois Institute of Technology: Drubka (1982), Shakib (1984), 

Corke (1987), Corke, Shakib and Nagib (1988). 

The observation that downstream large-scale vortex structures in a 
mixing layer, even turbulent ones, may modulate the primary instability 
near the lip came from Dimotakis and Brown (1976). Laufer and 
Monkewitz (1980) focused on the feedback due to pairing, but the 
complex coupling loop in Fi g 7 was deciphered at IIT. The reader is 
referred to Ho and Huerre (1984) for a rewarding description of the 
rich variety of other stability-related phenomena in initially laminar 
mixing layers. The IIT evidence consists of simultaneous measurements 
of unsteady pressure on the vertical flange above D in Fig 7 . velocity 
fluctuations sensed by a hot-wire movable along and across the layer 
and of various combinations of processing of both signals. Although we 
shall outline the developments in terms of single -wire data in the 
inset on the right of Fig 7 . the chicken-egg reconstruction of the 
system rested ultimately on two -point correlation, coherence and 
bicoherence functions of pressure and velocity. 

At the start of the layer in a sound- insulated low- disturbance 
2 

facility the (u') velocity spectra are broad-band, almost without 

peaks, indicating absence of spurious acoustic or vibrational 
disturbances and the weakness of the feedback. The rise of the first 
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two feedback- initiated frequency peaks is sketched in the inset of 
Fig 7 . Both the logarithmic slopes agree well with the linearized 
spatial amplification rates computed by Monkewitz and Huerre (1981) for 
the primary instability. The puzzling fact that the strongly amplified 
fundamental of frequency f is generally 15-25% slower than the 
theoretically strongest amplified frequency f was first mentioned to 

the author by A.K.M.F. Hussain ten years ago. Also initially 
surprising was the fact that its subharmonic at f/2 started at a level 
about an order of magnitude higher , see D in inset. The probable 
answer appears to be first that the nonlinear interaction of the 
combinations may extract energy more efficiently from the mean flow 
than the f^, f m /2 combination. Secondly the raltively high initial f/2 

level stems from the fact that vortex pairing is vigorous and generates 
higher pressure fluctuations than the more gradual vortex roll-up of 
the primary instability. The proximity of the latter to the lip is 
compensated for, since the pressure decay from a pairing multipole 
scales with wave length and therefore lower frequencies reach farther. 
The low Mach number equation for generation of the feedback pressure is 
shown at the top of Fig 7. The conversion into vorticity fluctuation 
is proportional to the spatial gradient of the pressure around the lip, 
not the pressure level, Morkovin and Paranjape (1971). Sharpness of 
the corner causes a large local increase in 3p/ds, thereby enhancing 
receptivity at separation. 

Up to the position A*, f and f/2 develop independently, propagating 
at different speeds close to those of linear theory. At A* a nonlinear 
threshold is reached where the subharmonic locks onto the speed of f 
and can thereby extract energy through the fundamental. This 
corresponds to the onset of the secondary parametric instability of 
Section 2.3: the subharmonic grows now exponentially at the linearized 

rate associated with the new base flow of rolled up vortices A and B of 
the fundamental. The subharmonic saturates at C, roughly at the 
position where the two strained vortices pass each other in x and 
manifest largest gradients du^/dx^ before viscous vorticity diffusion 
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cuts them down. (The products of the gradients in the source for the 
pressure fluctuations are summed on i and j . ) 

This double -frequency, globally unstable system was called weak 
because the fluctuation level at separation remains low, roughly 
hundred times smaller than that at the Nishioka-Sato separation on the 
cylinder, at Re-120, described in Section 3.2. Edge-tone systems also 
pump themselves up to higher energy levels since dipole impingement 
sources are more efficient than the quadrupole sources in the pairing 
process of Fig 7 . The strength of this compound mixing- layer 
oscillator is clearly limited by the saturation levels it can sustain. 
It can be disorganized by environmental disturbances and by artificial 
pressure fluctuations once these exceed sufficiently the initial levels 
of the fundamental and its subharmonic at D. 

(4.2^ Noise amplification in mixing layers 

The coexistance of independent linear development of f and f/2 in 
the range DA* illustrates the important fact that the mixing layer can 
amplify s imul t ane ous 1 v (albeit at different rates) frequencies from 
zero to about U/2 8 if imprinted at D. Only when the most prominent 
oscillation reaches sufficient nonlinear levels it may strengthen (or 
inhibit) other frequencies, such as the subharmonic at A*. This is 
effected through higher-order instabilities, often of parametric nature 
when frequency ratios are rational numbers, see Section 2.4. Further 
details on the richness of such forced nonlinear behavior in mixing 
layers are found in Ho and Huerre (1984) . 

Two consequences follow. First, low-frequency vorticity waves 
often permeate the shear layer, imprinted upon it at the origin by 
stream turbulence, irregular pressure fields, and feedback from 
large-scale downstream turbulent structures, Dimotakis and Brown 
(1976). At low amplitudes, their structure is little modified by the 
smaller-scale secondary instabilities, which they "modulate". After 
the layer thickens enough for them to acquire higher amplification 
rates they suddenly become active in a higher-order instability of 
their own. Sutton et al (1981) demonstrated the effect by vibrating 
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the separation point; the behavior persisted into the turbulent regime! 
It seems that a mixing layer, in contrast to the boundary layer, almost 
never forgets its long-wavelength vorticity structure. 

The second consequence is that extrinsic noise in a mixing layer 
will generate amplified, mildly filte red, noisy structures in 
accordance with the Ginsburg- Landau model of Deissler and Kaneko 
(1987). Thus, because its global instability is weak and limited to a 
short initial stretch, a mixing layer behaves quintessentially as a 
convectively unstable system, even when turbulent; see also Fiedler et 
al (1985). If a strange attractor existed for this x-dependent layer, 
its intrinsic irregularity would be contaminated by the extrinsic 
stochastic noise amplified by the layer. 

(4 , 3) Illustration of unsteady forcing of boundary layers 

In Section 2 we discussed unstable vorticity restructuring, 
assuming that vorticity nonhomogeneities were already internalized in 
the shear layer. In Section 4.1 and Fig 7 we described the manner in 
which an unsteady pressure gradient is converted into an unsteady 
internal vorticity wave at the origin of a mixing layer. This is a 
case of a strong very localized receptivity, proportional to dp/ds 
around the lip, with coefficients depending on the exact geometry and 
the forcing Strouhal number, f$/U. Next we examine a layer without any 
separation or two d imens i ona 1 leading edge . 

Figure 8 . courtesy of Kegelman (1982 thesis at Un. of Notre Dame) , 
displays visual evidence of the excitation of TS (Tollmien- Schlichting- 
Schubauer) , the viscosity- tuned waves on an ogive-cylinder model (see 
inset) by sound from a loudspeaker upstream of an open- throat wind 
tunnel at levels of 105-115dB. The smoke from an upstream source 
impinged symmetrically on the sharp nose and made the forced unstable 
waves visible some distance past the shoulder. Diagnostics with a 
hot-wire anemometer showed that the self -excited wave -pattern response 
to the forcing sound field began building up from the oscillatory 
Stokes layer -just upstream of the ogive -cylinder junction , where no 
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imprint on the smoke is yet visible. The range of excitation 
frequencies from 551 to 785 is impressive for an attached boundary 
layer. The spanwise inhomogeneous break up, featuring A vortex 
formations are now identified as H and occasional K secondary 
instabilities, Herbert (1988), as described in Section 2.3. 

Figure 8 demonstrates beyond doubt that harmonic forcing by 
unsteady irrotational fields generates early unsteady rotational 
waves in convectively unstable shear layers. Clearly, transition to 
turbulence cannot be predicted without the knowledge of the forcing 
field (frequency, amplitude, orientation), the dominant receptivity 
mechanism, and the threedimensional inhomogeneities of vorticity 
internalized in the boundary layer from other environmental 
disturbances. Even though the Reynolds number is the same in Fig 8 . 
the turbulent boundary- layer structures over the body appear to be 
governed by the above factors . Any more universal asymptotic features , 
if they could develop, evidently would require longer evolution. 

Of the three factors limiting the predictability of transition to 
turbulence, the receptivity mechanisms to harmonic pressure forcing 
such as that in Fig 8 have recently undergone considerable 
clarification. They will be discussed in conjunction with Fig 11 in 
Section 4.6. 

(4.4) On forcing of shear layers bv ouasisteadv fields; Goertler_ 
instability 

Forcing by steady vorticity fields and surface roughness has 

been largely ignored in the literature. And yet it is the primary 
source of the threedimensionality in secondary instabilities. 

Composite Figure 9 . based on wind-tunnel experiments of V. Kottke 
(1986) , illustrates the forcing of steady- in- the -mean Goertler vortices 
in a concave boundary layer (right side of Fig 2 ) by a primarily 
unsteady vorticity generator , namely the grid with mesh M and rods of 
diameter d in upper Fig 9 . Kottke (1980) and colleagues, e.g. Kottke 
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and Schmidt (1985) developed a chemical method for flow visualization 
in a wind tunnel and simultaneous photometric determination of local 
mass transfer rates. Grids and honeycombs bring forth complex 
threedimensional, multiply unstable and x-decaying wakes, Loehrke and 
Nagib (1972), (1976), which have been used for both generation as well 
as decay management of freestream turbulence. The curvature of 
Kottke's plate is relatively large so that in the parameter range of 
Fig 9 its boundary layer remains a buffeted laminar one , Section 2.6. 
However, the boundary layer is susceptible to centrifugal instability 
and responds to the mean- flow nonuniformities . which are invariably 
present (though seldom mentioned or documented) in flows with higher 
turbulence . 

The second row of Fig 9 demonstrates how the streamwise vortices 

form readily when the grid distance x g from the plate leading edge is 

small. As x increases the mean streamwise vorticity u> entering the 
s x 

concave boundary layer (or induced within it) subsides and so does the 

strength of the Goertler vortex response. In the third row of Fig 9 . 

the grid distance is held fixed, but the mesh M is varied, eliciting a 

wave number response in limited qualitative agreement with "parallel” 

linear theory for M> 7. The quantitative inference on the right of 

Fig 9 of the local and (spanwise) mean mass transfer coefficients 0 and 

in m/h supplement the visual information, 
z 

Clearly, the forcing by mean "environmental" » x distributions feeds 

the development of the centrifugal instability and predictions without 
better knowledge of the oncoming disturbances and the associated 
receptivity paths will prove fictional and wishful. For this flow, the 
critique _of_the formulation of the linear lnstabilitY_the_o_ry itself by 
Hall (1982), (1983) commpounds the predictive difficulties. 

Fortunately, recent and ongoing judicial application of the numerical 
marching technique of Hall by Kalburgi, Mangalam and Dagenhart (1988) 
and others restores considerable engineering utility to the results of 
the linear theory of Floryan and Saric (1984) . 

Amplitude resolution of measuring techniques for steady u> x fields 
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is orders of magnitude inferior to measurements in unsteady fields 
which are aided by the electronic amplification. Thus the comparison 
between linear theory and experimental results (many detected only at 
the nonlinear level) has always remained ambiguous. The net evidence 
points to strong imprinting by the particular distribution in each 

facility. In other words, the system is not robust with respect to 
disturbances, despite Huerre's 1987 conjecture that the system is 
absolutely unstable. 

In the experiments of Kottke in Fig 9 . only the effects of the 
primary instability are evident. Secondary unsteady instabilities, 
symmetric, sinuous and others have been studied experimentally only 
recently, and are burdened by the uncontrolled nonuniformities of the 
saturated primary field. For some results and guide to earlier 
literature see Swearingen and Blackwelder (1987). 

Figure 10 . middle, singles out two often neglected sources of 
quasisteady forcing of boundary layers: nonuniformities in leading 

edges (especially at supersonic speeds) and surface protuberances or 
dimples. Obviously the resulting local pressure gradients cause local 
cross-flows and At the author's suggestion Wilkinson placed a 

small intentional roughness element on the surface of a rotating disk 
and demonstrated in great detail that the well-known growing spiral 
vortices were driven by it, Wilkinson and Malik (1985). This forcing 
experiment carries special significance because Mack (1984) was able to 
reproduce to a remarkable degree the wave patterns and growth on the 
basis of linear theory, thus verifying this path of receptivity. 

Quasisteady forcing of instabilities is equally important in free 
shear layers: any disturbance entering the layer at the separation 

line where it is bom governs subsequent evolution of its 
threedimensional ity . Recent evocative visualization and discussion of 
the effects can be found in Lasheras and Choi (1988). 
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(4 .5) Environmental forcing in open-flow systems 

Most open- flow systems of technological interest are convectively 
unstable, Section 3 and Fig 6 . Therefore small environmental 
disturbances do energize and amplify the numerous competing intrinsic 
unstable modes in all primary and higher-order instabilities leading to 
turbulence. These external disturbances can be classified, Fig 10 . as: 
(A) irrotational unsteady pressure fields (sound, including near- field 
of pseudosound) which travel across streamlines, (B) threedimensional 
vorticity fields (usually decaying remnants of turbulence from 
upstream shear layers), convected with the mean motion of the fluid, 
while nonlinearly diffusing; (C) unsteady threedimensional 
entropy- temperature -density (ETD) spottiness of Kovasznay (1953) also 
convected "parabolically" from upstream mixing of shear layers of 
unequal temperatures; (D) steady- in- the-mean streamwise vorticity 

and ETD patterns, (present in all ground facilities but undetectable 
directly because of the exceedingly poor resolution of detectors; 

these fields force cross-flow instability and Goertler instability 
responses; (E) particulates and aerosols (including dust) in the fluid 
which individually and collectively cause flow nonhomogeneities along 
their nonequilibrated paths through the shear layer, e.g. Lauchle and 
Gurney (1984). In addition, (F) vibrational motions of any walls in 
the system and (G) nonhomogeneous or unsteady heating of such walls 
often provide significant seeding of instability modes (which may 
explain some anomalies in past observations). 

Because they are generally small, the disturbances are very hard to 
identify. Even with two independent instruments, the distinct 
irrotational velocity (pressure) fields and the rotational (laminar or 
turbulent) velocity fields have never been adequately resolved in any 
facility. The best calibrations of wind tunnels, such as that of 
Winter and Maskell (1980) leave hosts of unanswered questions as to the 
actual forcing fields. For practical purposes we don't know and are 
unlikely to know the input freestream distrubances into our unstable 
shear layers. In fact, they are not fixed and change during any test 
and from test to test. Philosophically, this reality must be 
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accomodated in any rational "prediction scheme" - perhaps by plausible 
bracketing estimates. 

The freestream fields, however, undergo straining and distortions 
by any solid boundaries, Fig 1 . Receptivity paths are to be sought with 
respect to the strained fields , which often are more effective. In 
fact, it is the straining of the sound field in Fig 8 . which enables 
the seeding of TS waves despite the disparity in wavelengths between 
the sound and vorticity waves of the same frequency. The straining 
provides new characteristic lengths commensurate with the TS 
wavelengths. Much remains to be done in the matter of straining of 
vorticity and sound fields as related to receptivity paths. At 
supersonic and hypersonic speeds the interaction of these fields with 
body shock waves creates additional unsolved problems. 

Various other issues and problems are listed in Fig 10 : they are 
essentially self explanatory . The overall perspective provided by Fig 
10 is sobering with respect to predictability of turbulent onset in 
open-flow systems. One implication points to the futility of judging 

the environmental disturbences by the single parameter \i /U^ at a fixed 

point, which is incorrectly identified as turbulence level. General 
correlations of transition based on this number have the scientific 

validity of schemes for choosing lottery numbers. When u / /U tfj exceeds 

0.5% or so, prospects of establishing reliable trends with various 
parameters are slim. Levels below 0.1% (including low-frequency 
contributions) are considered necessary for serious research on 
transition . In contrast experiments on controlled primary 
instabilities are considerably less sensitive to external disturbances. 

(4.6) Receptivities to u nsteady pressure, .fields 

The conversion of environmental unsteady pressure fields into 
unstable vorticity waves, such as displayed in Fi g 8 , has been subject 
to numerous experiments, with rather contradictory conclusions. The 
more significant experiments were summarized and critically reviewed in 
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Sections 1 and 5 of Nishioka and Morkovin (1986). The same paper 
describes the results of a more controlled experiment on acoustic 
irradiation of a Blasius boundary layer at unstable, neutral, and 
stable TS frequencies. As noted in the preceding section, the mismatch 
between the TS and acoustic wavelengths (or equivalently between the TS 
and sound propagation speeds) tends to make this receptivity rather 
poor. 

Receptivity improves as additional effective characteristic 
lengths are introduced into the combined system of the forcing, 
diffracted acoustic field and the boundary layer. An infinite 
Poiseuille-duct flow has zero linear receptivity to a uniform harmonic 
pressure gradient in the duct. Both flow solutions, being uniform in 
x, simply superpose without interaction. Variations of the amplitude 
of the unsteady pressure gradient with x, as in Fig 8 . and variation of 
the boundary layer thickness with x provide two specific receptivity 
paths; see Nishioka and Morkovin, Sections 2 for the physics of the 
interaction and Section 5 for other receptivity paths. Because of the 
lack of controlled experiments, the question of which paths are most 
efficient for a given geometry remains open. 

Depending on geometry the seeding of TS waves may take place over a 
short segment of the boundary layer, shorter than one TS wavelength. 
Such localized receptivity paths lend themselves to efficient analysis 
based on the triple-deck approach; see review by Goldstein (1986). On 
the other hand, extended distributed contributions to the TS vorticity 
strength do occur frequently over bodies with continuous curvature. 

The analysis is then complicated by positive and negative phase 
interference effects. 

The left side of Fig 11 outlines the formal solution for the 
forcing, i.e. for the nonhomo geneous differential equations system. In 
principle, it depends on the completeness of the eigenfunctions for the 
system, such as proved for the Orr-Sommerfeld equation by Salwen and 
Grosch (1981). It was carried out explicitly by Bechert (1988) for the 
simpler inviscid case of a vortex sheet downstream of a semi- infinite 
plate and buttressed by the low-frequency experiments of Bechert and 
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Stahl (1988). 

The familiar resonant vibrations of a wine glass in response to the 
excitation of a singer's voice will illustrate both the similarities 
and differences with the forcing behavior in an unstable boundary 
layer. The forcing fields evoke particular solutions of the 
nonhomo gene ous linear differential equations governing the forced 
system. To these we can add any number of solutions of the homogeneous 
unforced system with arbitrary coefficients, the normal modes or 
eigensolutions . The requirement that the superposed particular and 
normal-mode solutions satisfy the complete boundary conditions over the 
boundaries of the wine-glass or the shear layer generally fixes the 
coefficients of the eigensolution spectra and thus selects 
quantitatively the response to the excitation. 

In the case of the wine glass the vibrational eigensolutions are 
all frictionally damped and the response consists primarily of the 
particular solution; it is essentially proportional to the sustained 
amplitude of the singer's tone. Once the tone ceases, the response 
rapidly decays. In the rare cases when the resonant response amplitude 
exceeds a nonlinear threshold, the properties of the physical system 
are no longer governed by the idealized linear differential equations 
and the glass breaks . 

The corresponding Nishioka-Morkovin experiments with a very weak 
sound source in proximity of a Blasius boundary layer demonstrated that 
when the TS waves at the exciting frequency were damped according to 
the linear Orr- Sommer f eld equations, a behavior ensued, similar to that 
of the wine glass. At each x, z, the response consisted of a 
finite -wavelength oscillatory viscous Stokes field in y, (the 
particular solution) plus small contributions from the damped 
eigensolutions. However, when the forcing frequency was in the 
spectral domain of amplified normal modes, the cumulative response grew 
in x very much faster (!) than allowed by the theoretical amplification 
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rate of the unforced normal modes. Once locally excited, these TS 

waves propogate downstream and continue to grow even if the periodic 
excitation ceases - in sharp contrast to the case of the wine glass. 
This is of course the behavior expected in a convectively unstable 
layer, inferable from Fig 5 . As in the case of excitation by a 
vibrating ribbon, the propagation field becomes that of the amplified 
homogeneous normal mode. The key difference is that its amplitude is 
not arbitrary (as in linear stability theory) but set linearly by the 
amplitude and geometry of the forcing sound field, i.e. by the net 
receptivity for the configuration. 

When the unsteady pressure field is not periodic, the linear 
mathematics becomes more involved, but the some principles govern the 
seeding of the amplified disturbances which are convected downstream 
within an unsteady laminar "spot". Kendall (1987) induced a traveling 
pressure wave as a cylinder approached and receded from a flat-plate 
boundary layer. The resulting complex amplifying pattern was 
documented through simultaneous measurements with multiple pressure 
sensors at the wall and hot-wire anemometry in the layer. The findings 
are consistent with the concepts described in this section. 

(4.7) Unknown receptivity paths 

When the external disturbances are large, additional nonlinear 
receptivity paths appear. They remain uncharted except for a few cases 
of high-level sound excitation, Ginevskii et al (1978). The bypass 
roads to turbulence, discussed in Section 2.6, by definition fall in 
the category of unknown receptivities. The right side of Fig 11 . lists 
the more important cases and conjectures that local and/or transient 
inflectional shear-layer profiles may be responsible for the rapidity 
of the turbulence onset. Gill (1965) calculated that such an effect 
could be responsible for the transition in the linearly stable plane 
Couette flows and pipe flows. 

However, unknown linear receptivity paths provide the conceptually 


more important challenges . Although Fig 9 demonstrated forcing of 
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Goertler instability by quasisteady distributions, the corresponding 

linear receptivity mechanisms remain unknown. Similarly the forcing 
paths of the steady cross -flow instability as well as of the 
threedimensional traveling waves in threedimensional boundary layers 
have not been discussed in the literature. 

Perhaps the most vexing and technologically important is the puzzle 
of receptivity to low freestream turbulence . Both experimentally and 
theoretically this is the most difficult linear receptivity problem. 

As noted on the left of Fig 11 . any analysis performed using averaged 
quantities would solve the wrong problem. Almost of necessity, 
experimenters must use space- time correlations and coherences between 
multiple sensors. The difficulties far exceed those of known traveling 
wave fields of Kendall (1987) at the end of Section 4.6. 

Hot-wire evidence: Bennett (1953), Arnal and Juillen (1978), 

Klebanoff (unpublished, 1964) and intriguing visualization in adverse 
pressure gradients by Gates (1980) imply that the receptivity path for 
low- intensity turbulence includes TS wave-packet formations, as the 
convective nature of the instability would suggest. Kendall's more 
controlled experiments (1984, 1985) with more extensive up-to-date 
instrumentation do indeed disclose energy in the TS bands (with 
decreasing frequency of the broadband peaks as the layer grows in x) . 
However, serious questions remain whether this path represents a weak 
ineffective sideshow while some threedimensional instability governs 
the dominant road to turbulence. 

For further discussion of two possible ways of internalization of 
the external vorticity disturbances see p. 256 of Nishioka and Morkovin 
(1986) based on the implications of the cited research of 
Arnal -Juillen, Klebanoff, and Kendall. It is quite possible that the 
receptivity paths may depend on the specifics of the oncoming 
low- intensity, decaying turbulent wave packets in any given 
realization and on the straining fields associated with the flow 
facility and the front of the body including those at sharp leading 
edges . 
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Finally, a curious hybrid case is noted at the bottom of Fig 11 . 
Unexpectedly, moderate distributed roughness which might disorganize TS 
waves, accelerated (!) the TS wave formation and secondary 
instabilities, Corke, Bar-Sever and Morkovin (1986). These authors and 
especially Kendall (1981) showed that the boundary layer just beyond 
the roughness peaks follows displaced Blasius profiles, i.e. is 
governed by molecular rather than eddy transport. The plausible 
inference from the total evidence suggests that the low- inertia flow 
regions below the roughness peaks modifies substantially the 
receptivity to external disturbances from that of a smooth wall. When 
the TS wavelengths exceed 15-20 times the characteristic roughness 
scale, the crucial y and x phasing of TS waves evidently remains 
unimpaired. These inferences should be verified by controlled 
experiments designed for the specific testing of the mechanisms. 

The preceding roughness role in transition is distinct from the 
larger roughness role covered on the right side of Fig 11 . and was not 
previously suspected. The third technologically all-important 
roughness role on the roads to turbulence is that of very small 
distributed roughness . This is probably not a receptivity problem: at 

such low Reynolds numbers, the roughness generates steady vorticitv . 
especially Any unsteadiness must be associated with an instability 

forced by a "cooperating” disturbance agent, such as freestream 
turbulence. Conceptually, the small -roughness effect remains an 
irritating bypass , almost an insult to our engineering soul. 
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5. CHAOS AND TURBULENCE IN OPEN -FLOW SHEAR LAYERS 

(5.1) Instabilities, forcing dis turbances, receptivities and chaos 
prediction 

In Section 2, we outlined the observed development of sequences of 
instabilities and onset of turbulence in x-dependent open-flow systems. 
In Section 3, we established that except in the infrequent cases of 
global instability, environmental disturbances (not the intrinsic 
properties of the initial shear layer) control the instability 
sequences on the road to turbulence, and dominate the irregular 
character of vortical structures. Thus in this vast class of 
technologically important flows the concepts of chaos and strange 
attractors offer little if any help with the design task of predicting 
transition. 

Since, as explained in Section 4.5 in conjunction with Fig 10 we 
cannot expect to have adequate information about the environmental 
disturbances - the input into the instabilities - we must establish 
rational bracketing estimates in spite of these uncertainties and a 
matching philosophy of risk in design; see Morkovin (1978) for further 
discussion. For the purposes of the present review, only the 
rationally solvable aspect of the problem of receptivity paths for the 
external disturbances were taken up in Section 4. That completes the 
account of the main recent conceptual and physical advances in 
understanding transition to turbulence In open- flow systems. 

In the Appendix, we provided the essential descriptive background 
for the conceptually and philosophically important notion of chaotic 
behavior. We can now return to the questions posed in Section 2.5 
concerning the relation of chaos, transition to turbulence and 
turbulence in open- flow systems. In Sections 5. 2 -5. 4 we shall 
reexamine three flows (one of them globally unstable) to set the stage 
- for more questions. Currently the uncertainties and irregularities 
of the subject in open- flow systems lead naturally to questions and few 
answers if any. However, with the recent clarification of concepts of 
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convective and global instabilities, the questions should be sharpened. 
In fact, an experimenter should well consider the various questions in 
Section 5. 2 -5. 5 and Fie 13 and 14 before embarking upon the 
time-consuming search for chaotic characteristics in open- flow shear 
layers . 

(5.2) The case of a chaotic laminar slender wake 

The beautiful visualization of Fig 12 of S tuber and Gharib (1988) 
display for us four modes of behavior of a slender wake downstream of a 
63A008 airfoil under various intentional excitations by local wall 
heating in water. These should help to make the concepts in the 
Appendix and Fig A-l more concrete - and raise some more general 
questions stimulated by the behavior of this specific fluid oscillator. 

The slender wake is convectivelv unstable so that it responds to a 
band of frequencies which widens with excitation level (though the band 
is much narrower than that of mixing layers). Subject to the 
indigenous environmental disturbances the spectrum at the primary probe 
position 3 chord lengths downstream shows a broad-band hump between 2.2 
and 5.8 Hz, with most energy between 4 and 4.6 Hz. This is in sharp 
contrast to the globally unstable bluff wakes and mechanically rigid 
oscillators, where the oscillations are very regular. Nevertheless the 
investigators identify 4.075 Hz as the "natural" frequency f and study 
the properties of this nonlinear Navier- Stokes system with rather 
strong forcing, as seen in Figs 12a-c . 

The forcing at f^ in Fig 12a makes the oscillator more regular: 
the excitation overrides the competing frequencies in the convective 
instability and reduces the background broadband spectrum by factors 
roughly from 10 to 50. The quasiperiodic f^ and f^ forcing of Fig 12b 
should place the system on the torus of the 3 phase space in Fig A-lc . 
although spectra in Fig 7 of S tuber and Gharib (1988) indicate that the 
oscillations are more complex. Evidently the convective instability 
elicits more degrees of freedom in the open- flow environment. 
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In accordance with Section A- 3, forcing by the three frequencies in 
Fi g 12c should make the wake chaotic (as the spectra seem to confirm) . 
However, according to S tuber and Gharib (1988) examination of more 
extensive visualization records disclosed "some areas of strong 
mixing". Since in Fig 12c no local fuzziness nor obliteration of 
streaklines (Section 2.5) is evident, the original label of "laminar" 
was retained for this specific realization. In principle, 
twodimensional laminar chaos should be possible and may presumably 
correspond to so-called twodimensional turbulences (Section A-4) . If 
at this speed the local turbulent mixing is sporadically present, a 
reduction in Reynolds number would probably produce a true laminar 
though not necessarily twodimensional chaos. 

Sporadic local onset of turbulence is common in bluff wakes. In 
fact, its march upstream towards the body with increasing Reynolds 
number modifies the wake characteristics significantly, Fig 13 . top. 
Such localized onsets of turbulence in a sea of laminarity bring out a 
basic question: how does one usefully characterize the state of a 
shear layer in the theory of dynamical systems when it displays 
laminarity and turbulence side bv side ? See also item (2) on the right 
of Fig 1* . 

The visualization in Fig 12d was carried out by K. S tuber at the 
author's request. The turbulent structures are "natural" in the sense 
that no oscillatory forcing was applied in this case. The boundary 
layer was turbulent at the trailing edge of the airfoil. The fuzziness 
of the lines testifies to small-scale turbulent diffusion which is not 
present in the other visualizations. The intent was to suggest that 
turbulence and chaos are not synonymous in open- flow systems. 
Paradoxically, the same conclusion was reached recently by Heslot, 
Castaing, and Libchaber (1987) in their experimental study of a 
Benar d - Ray 1 e i gh cell - a closed-flow system ! Further discussion of 
this issue will be taken up at the end of Section 5.5. 

The paper of Stuber and Gharib also illustrates the nature and the 
approximate validity of temporal phase -space reconstruction and 
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Poincare" sections for the cases of Fig 12a-c . However, the 
determination of Lvaounov exponents (Section A- 3) in their Fig 17 shows 
them to be positive for the non- chaotic states of Fig 12a and b . This 
is a significant conceptual discrepancy which they "attribute to noise 
inherent in an open system". Clearly the value of the concepts of 
chaos deteriorates unless rigorous procedures can be established for 
&Q X . s e para ti ng, £h-g.. l n£XlraiS- £ ,h§.9..Ltc -featu re s from extrinsical driven 
random features in the flow response; see item (3) on left of Fig 14 . 

All the measurements were taken at a single point or single x 
station. How well can such measurements characterize the properties of 
the whole system, especially in an x-developing shear layer? (Item ( 2 ) 
on left of Fig 1 . ) Admittedly, Deissler's reasoned call for 
measurements of Lyapunov's exponents in a convective frame of 
reference . Section A-4 . is difficult to carry out experimentally. 
However, the "optimal space-time correlations" of Favre et al (1953) 
substantiate Deissler's criticism. Since they take into account the 
convection of large structures they may lead to improved procedures: 
see Question in Fig 14 and Section 5.5. 


The preceding questions sparked by the Stuber-Gharib experiment do 
not imply any criticism of the experiments. After all their intent was 
to "test the waters" for the applicability of the new concepts and 
procedures in a relatively simple and controllable open- flow system. 
They accomplished their purpose well. The author is grateful for their 
help with Fig 12 and for an early manuscript, which helped to broaden 
the range of the planned general questions. 

(5.3) Road to turbulence in a globally unstable flow 

In view of the central role of convective instability in forcing 
and controling the sequence of instabilities on the road to turbulence, 
it is instructive to examine the major features of a similar evolution 
as Reynolds number increases in an unidealized globally unstable 
open- flow system, namely that downstream of circular cylinders, upper 
part of Fig 13 . Both the globally unstable near wake as well as the 
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local incipient turbulent spots at some Reynolds numbers. Therefore, 
the two questions at the top of Fig 13 are in order in the context of 
Section 5. 

In a globally unstable near wake vortical structures are still 
convected away without the chaos -generating iterative recycling present 
in closed- flow systems. However. part of their unsteady signature is 
conveyed upstream and thereby helps to create a partially closed system 
with respect to temporal disturbances . Thus open- flow globally 
unstable systems have one important feature in common with closed- flow 
systems. The applicability of chaos concepts may therefore be more 
likely in globally unstable cases than in the convectively unstable 
ones . 

Sections 3.2 and 3.4 established the basis for the global 
instability of near wakes and described the temporal evolution to the 
final saturated x- resonant flow field. Measurements by Nishioka and 
Sato (1978) and Unal and Rockwell (1988) in this final field disclose a 
spatial exponential growth in x of the harmonic oscillations u' . The 
strict global -instability regime extends only between (aspect-ratio 
dependent) Reynolds numbers of 45 and 160, with two minor puzzles, 
still unresolved. The "clockwork oscillatory regularity" of Roshko 
(1955) is interrupted by small discontinuities near Re of 90-100, 
apparently associated with a change of shedding pattern and possible 
threedimensionality , see discordant accounts of Tritton (1959), Gaster 
(1971) Berger and Wille (1972), Gerrard (1978) and Friehe (1980). 
Furthermore, in some low-Re experiments, the phase of the vorticity 
varies linearly along the span yielding swept-back vortices. 

Smaller- scale spanwise nonuniformities appear near Re of 200. 
Gerrard' s (1978) visualizations suggest the presence of vorticity 
feedback to the separation line - as postulated in absolute instability 
- for this spanwise threedimensionality, Fi g 13 . At about the same 
Reynolds number local regions of high diffusivity (incipient turbulence 
spots of Section 2.5) appear at first quite far downstream in the 
convective- instability domain. Also at about the same Re, a startling 
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change in the nonlinear upstream influence takes place: a reversal in 

the Re -trend of base pressure and oscillatory lift coefficient on the 
cylinder (both of technological importance) ; see introduction in Unal 
and Rockwell (1988). These trends reverse again for Re > 2000-10,000 
(depending on aspect ratio) . 

Despite various conjectures, Gerrard (1978), Unal and Rockwell 
(1988), no convincing explanation has yet emerged for either reversal . 

A correlating scale seems to be provided by the formation length of 
Gerrard, 1^ sketched in Fig 13 . This length corresponds to the 
crossing of the centerline by the rolling-up vorticity sheet when the 
first vortex is shed. The system, which now includes threedimensional 
vorticity, evidently remains globally unstable. The cause of the 
reversals must therefore be sought in the changes of the vorticity 
distributions and (possibly competing) instability modes in the 
immediate near wake as Re increases. Symmetric Gerrard-Bloor 
instability (which at higher Re is associated with the thin individual 
separating layers rather than with D) and the aforementioned 
smaller-scale spanwise threedimensional instability could well 
interfere with the basic Karman-Benard mode. The simple global 
instability model probably will not be able to account for the subtle 
changes in these vorticity fields and their multiple instabilities as 
Re increases. Nor can one expect chaos theory to unravel the puzzle in 
this highly x- dependent shear layer. 

With growth in Re comes also the upstream march of the first onset 

of turbulence Tu, Fig 13 . We could expect changes in behavior when Tu 

first reaches the boundary between global and convective instability or 

the fluid in the first shedding vortex. These conditions probably 

modify but do not destroy the global instability of the main 

oscillator, but the author is unaware of any studies of such effects. 

The Gerrard-Bloor vortices increase in frequency and become smaller. 

At higher Re many of them crowd into the shear -layer segment which 

rolls up into the first large vortex at 1 . The first bursts of 

r 

turbulence then apparently become associated with higher instabilities 
of the strained unsteady Gerrard-Bloor vortices. They ultimately reach 
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the proximity of the oscillating separation line on the cylinder. 

It has been recognized only in the Fifties that the subsequent 
classical rapid decrease in drag takes place through the formation of 
laminar bubbles when the oscillating free turbulence shear layer 
reattaches to the cylinder. This terminates any global instability 
features in the near wake of circular cylinders. The reattachment 
process is random in space and time and very sensitive to freestream 
disturbances (especially in because the instability in the 
individual separated layers is convective . The entrainment effects 
within the bubbles are commonly uneven with consequent asymmetry of the 
mean flow. The rest of the listings in Fiq 13 for still higher 
Reynolds numbers are also technologically important. However they 
relate now to convectively unstable boundary layers under high 
acceleration and are included here merely for the sake of completeness. 

It seems clear that despite the important feature of global 
instability, the prospects of clarifying any of the significant details 
of these high-Re phenomena in an open- system with strong x- dependence 
through concepts of chaos are virtually nil. To be fair no such claims 
have been made for this specific applicability, as they have for the 
power of modern computers. It seems safe to repeat after 24 years, 
Morkovin (1964), the opinion that these too are unlikely to handle the 
important threedimensional and locally turbulent features of the flows 
described in Fig 13 for quite some time. 

It is worth remarking that much of the complexity of these flows is 
associated with the interplay between the unsteady, variable, often 
spanwise nonuniform, separation line and the nonlinear pressure 
gradients . In bluff flat -based wakes such as that of Hannemann and 
Oertel (1988), where the separation line is fixed at the straight 
shoulders global instability persists when the separating boundary 
layer becomes turbulent! If its thickness is small with respect to the 
base height, the high-frequency small-scale turbulence at separation 
has remarkably little influence on the wake oscillator including its 
frequency . 
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15 J&1 RoM-to turbulence in a boundary layer;, e N method 


The cases examined in Sections 5.2 and 5.3 had inflectional 
instability. In contrast the primary instability in the boundary- layer 
evolution depicted in the lower half of Fi g 13 is viscosity induced. 
Nevertheless, it is convectively unstable as was the slender wake of 
Section 5.2 so that many comments and questions need not be repeated. 
Since excitation by unsteady pressure gradients was covered in Section 
4.6, the response selected here for examination is that to low 
freestream turbulence. 

Despite Kendall's misgivings cited in Section 4.7 we assume that 
coherent freestream structures, strained and accelerated by the body 
pressure field enter the boundary layer or pass near it and induce 
laminar TS spots as mandated by satisfaction of all the boundary 
conditions, left side of Fig 11 . In accordance with the character of 
convective instability Fi g 5a . these spots grow while propagating 
downstream; they merge and superpose. A hot wire would measure 
modulated, slightly shifting frequencies of an x-amplifying u 
fluctuation. Since 2D waves are generally more amplified than skew 
waves, the averaged response would appear close to a 
quasi twodimensional growth in x. 

The sequence of 3D secondary instability, tertiary instability, and 
intermittent onset of incipient turbulent spots between Re ^ and 

Re . . sketched in Fie 13 summarizes the material discussed 

tr - end 

throughout Section 2. The label under the primary instability 
emphasizes the fact that the amplification rate of TS growth is slow 
and that the sketch distorts the x scales to accommodate the 
explanations. In reality the x-range between Re cr amd the beginning 
of the secondary instability vastly exceeds the combined x-range of the 
secondary and tertiary instabilities. 

A.M.O . Smith (1956) and Van Ingen (1956) suggested independently 
that this evident primary control of the transition process by the 
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involving TS instability. With the blossoming of computer technology 
this technique based on the ratio of linearized amplification between 

Re cr (f) for any frequency f and any station x, e v , has been widely 
applied and extended, e.g. Bushnell, Malik and Harvey (1988). The 
following comments on judicious use of the technique are appropriate in 
the spirit of Section 5. 

The crucial choice of N (say from 5 to 12 in free flight) clearly 
depends on estimated environmental disturbances for all convectively 
unstable shear layers. Philosophically, the risk to the design 
okje^tiYea and the possibility of bypass transition , fi g 15, in 
uncharted flow fields (including wall roughness generation in superhot 
hypersonic environments) should be accounted for. The very term 
"prediction" instead of "estimation" indicates a tendency to an 
optimistic overbelief by engineers: we tend to endow the results of 

our simplified computer codes for idealized flow fields with an 
unwarranted aura of reality. 

The basis for inferring the quality of "controlling slow- rate 
amplification" for primary instabilities other than the low- speed TS 
instability needs buttressing. The reliability of computer codes for 
calculating accurately enough subtle boundary layer profiles and the 
corresponding Re^ r (f) and amplification rates, integrated over the path 
of the laminar spots (e.g. for 3D boundary layers with wall curvature, 
hypersonic boundary layers, etc.) is largely untested. The critical 
"semiempiral" choice of N for any class of flows is a sub j ective 
process of judgments across several disciplines. It involves judgments 
across the Mach number and Reynolds number range on specific linear 
instability theories, on specific computer codes, on the onset of 
transition in the few experiments with sufficient detail information, 
on the quality of environmental disturbances in these experiments, etc. 

Such judgments have come largely from a single group at NASA 
Langley Research Center. As a friend and an old industrial hand at 
transition estimates, the author is supportive of the objectives and 
efforts of the group. Nevertheless, it is appropriate to quote here 
the modified Guideline No 4 of the U.S. Transition Study Group, 
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Reshotko (1976). "Whenever possible, tests and numerical computations 
should involve more than one facility. Tests and computations should 
have ranges of overlapping parameters, and whenever possible, 
experiments should have redundancy in transition measurements and 
transition-estimating personnel . w The issues here transcend the narrow 
professional questions about which author is more correct in what 
paper. In view of the magnitude of private and public funding 
involved, these are national issues. 

Figure 13 . bottom also cautions us about the so-called unit-Re 
effect. Whenever frees tream velocity or stagnation pressure is changed 
transition on the model changes not only in x but also in the 
diensionless parameter Re^! This pernicious nonconstancy of R« tr is 

usually conveyed in plots of Re vs Re per unit length, Re ; hence the 

tr JL* 

misleading name. The diverse variations merely indicate that for any 
model in any facility there are more than one characteristic length or 
velocity which influence palpably the transition process. There is no 
single unit-Re effect ! Notorious contributors to these variations come 
from leading-edge geometries, roughness and freestream disturbances, 
including sound radiated from turbulent sidewalls in supersonic wind 
tunnels. The unit-Re effects make it difficult to compare transition 
results between facilities and free flight and complicate the life of 
transition estimators. 

The comment on the right of the sketch raises an untackled issue 
for research in turbulent boundary layers : to what extent do the 

properties studied in a given realization reflect generic intrinsic 
turbulent- layer characteristics, uncontaminated by the nonuniformities 
in the preceding laminar layer including the extensive transition 
region? The question is relevant whether the transition occurs 
"naturally” or is artificially induced. 
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(5.5) On experiments on chaos in open-flow systems 

The account of observed real-life routes to turbulence in the three 
classes of flows in Sections 5. 2 -5. 4 prepared us for further discussion 
of the questions posed in Section 2.5. The central questions concern 
the manner in which the successive phenomena in the sketches of Fig 13 . 
could be related to dynamical states, including chaos: the definition 

of the state of the system (see Section 5.2) and the existence of a 
strange attractor associated with each segment of a shear layer which 
grows with x. The relative success of the linear instability theories 
rests on the fact that linear instability characteristics , especially 
the amplification rate, are almost-local properties of each segment of 
the shear layer. We can define the functional dependence a^(x) and 
use it in our computations. Can a strange attractor be related to each 
segment of a chaotic or turbulent layer growing in x by a specific 
rigorous procedure ? If not, repetition of techniques taken over from 
experiments on closed-flow systems, would be illusory with superficial, 
if any meaning. 

The notion that a major difference between open and closed systems 
lies in the fact that vorticity structures in the former are convected 
away and not recycled as in the latter was brought up in Sections A- 4 
and 5.2; see also the contrast between statements (a) vis-a-vis (b) and 
(c) in Fig 14 . Deissler (1985) in effect states that the divergence 
properties characterized by Lyapunov* s exponents are not 
"almost-local" and that a convective framework is needed to elicit them 
in his "simple" equation (i.e. with coefficients not explicitly 
dependent on x) . Are other properties . such as the Strange Attractor 
S.A., itself and its dimensions also attached to the traveling 

vo&ti£ % 

If that were so for Deissler' s model, determination of these 
properties would become questionable in all open- flows where the shear 
layers changes in x. The asymptotic target would be constantly 
changing so that for practical purposes the properties would be unde- 
f inable : see contrast in statements (1) on both sides of Fig 14 . 
Modernized spatio-temporal multiple-probe techniques of Favre and 
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coworkers, at bottom of Fig 14 . might possibly characterize the 
intrinsic properties in constant-scale turbulent flows such as those in 
nines . That would appear to be a meaningful but difficult experiment 
which should clarify the convective aspects of any chaotic properties. 
The studies of Wygnanski and Champaign (1973) of turbulent puffs and 
slugs provide support for this notion. 

Even pipe flows, Fig lb . contain segments of x-dependent internal 
boundary layer growth which make them at least locally convectivelv 
unstable . As such, pipe flows may be imprinted by extrinsic noise 
before they become fully turbulent. As in the case of the boundary 
layer at the end of Section 5.4, the question arises whether the fully 
turbulent structures in nines have anv memory of their noise-controlled 
past or how rapid the memory loss may be. This is related to the 
question of separation of intrinsic chaotic features from externally 
driven random features; see items (3) in Fig 14 and discussion of noise 
at the end of Section A-4. The turbulent pipe -flow experiment in the 
preceding paragraph could then be usefully followed up with studies of 
effects of diverse intentional disturbances which trigger different 
upstream transition patterns. 

Figure 14 summarizes the various issues which distinguish between 
concepts and their experimental verification in closed- flow and 
open- flow systems. The preceding discussion justifies the label 
"controversial" on top right of Fig 14 . and the skepticism with respect 
to single-point measurements, especially in open-flow systems. What 
substantive information can be inferred from a time trace at a single 
point of a growing turbulent shear layer? Only multiple -point 
multiple-probe information can help to clarify the modifications 
brought about by the uniterated convection of vortical structures. 

As noted in Section 5.2, the one experiment in the literature 
reporting two-point measurements is for a closed (!) Benard- Rayleigh 
cell, heated from below (in which buoancy- driven vorticity recirculates 
in place). The results created considerable commotion: Heslot, 

Castaing, and Libchaber (1987) confirmed that chaos and turbulence are 
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not synonymous. As the Rayleigh number increased they observed "the 
onset of the oscillatory instability followed by the now well-known 
routes to chaos", specifically through a period- doubling cascade 
proposed theoretically by Feigenbaum. The correlation dimension of the 
temporally chaotic state varied from 2 ± 0.1 to 4, before the coherence 
between the bolometer at half height of the cell and the bolometer near 
the bottom wall started decreasing from nearly unity for all 
significant frequencies to a negligible value. The loss of coherence 
and the disappearance of the identifiable chaotic state was attributed 
to the formation of a boundary layer at the wall. 

As the Raleigh number was pushed higher and higher they observed 
two extended regimes of "soft" and "hard" turbulence. The first, with 
very little coherence between the two bolometers, is interpreted in 
terms of unsteady rolls with locally detached laminar boundary layers 
and random phasing, a state called phase turbulence in some quarters. 
The hard turbulent state is interpreted in terms of a turbulent 
boundary layer at the base with abrupt detachment of thermal flumes 
which reestablish communication and coherence between the two 
bolometers (one above the other). Interestingly, the interpretations 
are supported by dimensional arguments from en gineering turbulent 
boundary layer theory , besides spectra, timetraces, probability 
density, etc. In other words, even in closed systems, we ultimately 
have boundary- layer formation, detached shear layers, etc. and need to 
speak of them as distinct entities, with turbulence on their local 
scales, as in Section 2.5. 

In the authors' language: "we have in some ways brought together 

the two approaches to turbulence. The dynamical system approach is a 
very good analysis of the onset of (temporal) disorder, as one slowly 
increases the Rayleigh number... It is followed by a state where s pace 
disorder sets in. In this regime, the development of a boundary 
layer, laminar at first and later turbulent, leads to two distinct 
turbulent states ..." 
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Chaos concepts are evidently of little help after the spatial 
disorder commences. Can the chaos concepts be at least conceptually 
fruitful in o pen* flow shear layers , where we begin with already formed 
boundary layers ? How much is their significance weakened by the 
downstream convection of vortical structures absent in the 
Benard-Rayleigh experiments. And by the growth in the layer thickness? 
Some of the suggestions in this section concerning experiments may help 
to provide at least partial answers. 
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6. CONCLUSIONS 

Roads to turbulence in open- flow systems can be interpreted 
fruitfully as sequences of often competing instabilities corresponding 
to primary and higher-order restructuring of vorticity distributions; 
see summary Figures 2, 3. 13 and 15 . The recently clarified concepts 
of convective instability and global instability unify previously 
fragmentary observations of fluctuation behavior, Section 3. In 
particular, a large majority of technologically important open- flow 
shear layers, which belong to the convectively unstable families, are 
unquestionably driven by ever-present small environmental 
disturbances. Consequently in these flows, concepts of chaos and 
strange attractors (introduced in the Appendix) offer no help with 
prediction of transition to turbulence. 

To help understand forcing by external disturbances the 
disturbances are first classified into families with different 
propagation characteristics, Fig 10 . These disturbances are first 
strained or diffracted by the wall geometry of tunnels and models. 
Depending on their propagation characteristics and the straining near 
and in the shear layer, these external disturbances are converted into 
internal unstable modes of Tollmien-Schlichting-Schubauer , Goertler, 
and cross -flow primary vorticity restructuring and into more complex 
modes of secondary and higher instabilities. The process of 
conversion was identified by the author in 1969 and called 
receptivity. Several types of receptivity paths are demonstrated by 
experiments selected from recent literature, Figs. 8___and 9 and 
Sections 4. 3-4. 6. When the external disturbances are small, the 
forcing corresponds to a nonhomogeneous linear problem, which is 
solvable in principle, Fig. 11 . The mathematical and physical basis 
of several receptivities is outlined. 

Receptivity analyses should alert us to the most dangerous 
receptivity paths for any given configuration. They will not help 
with direct estimates of transition because the disturbance flow 
fields will not be known a priori in wind tunnels or in flight. 

Because of the large total amplification of the disturbances, the 
initial forcing fields are too weak for adequate diagnostics and 
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measurements of the distinct families of disturbances of Fig. 10 . 

Since the input into the various instabilities will remain uncertain 
we have to resort to engineering estimates and bracketing with the 
powerful but tricky aid of computers. The philosophy of approach to 

N 

such estimates, in particular, the e method, is examined in Section 

5.4 In Section 5 we return to the questions posed in Section 2.5: 
To what extent are the classical turbulent states in open- flow systems 
equivalent to chaotic states? What deeper insights can the theories 
of chaos bring to our dealing with turbulence in open- flow systems? 

The attempts at answers raise new questions about the effects of 
vorticity convection and of x- dependence of the shear layers on the 
nature and existence of strange attractors in open flow systems - see 
Fig. 14 . The roads to turbulence in three real-life shear layers, 

(one with the partially closed system of disturbances rather than of 
mass and vorticity) are examined for clues in Sections 5. 2- 5. 4. 

Doubts arise about the meaning of single-point measurements taken over 
from experiments in closed-flow systems. Modernized multiple- 
probe optimal space -time techniques of Favre and co-workers are 
suggested to help to answer the vorticity-convection and 
x-dependence issues . 

The new experiments of Stuber and Gharib (1988) in slender wakes 
and Heslot, Castaing, and Libchaber (1987) in a closed Be^ard-Rayleigh 
cell heated from below suggest that chaos and turbulence are not 
synonymous. The latter, closed- flow experiment (the only substantial 
one yet reported using measurements at two points) indicates that 
chaos concepts are of little help once space disorder as contrasted to 
temporal chaotic disorder sets in. In this closed-flow, space 
disorder appears when the controlling parameter is high enough for the 
formation of boundary layers. In open- flow systems we have boundary 
layers and separated shear layers as an initial given! Thus current 
evidence points to negative answers to the two questions about chaos 
and turbulence . 

In summary , the concepts of convective and global instabilities 
help us to unify diverse observations on instabilities and transition 
in open- flow systems and confirm the role of environmental 
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disturbances. The roads to turbulence are physically and 
theoretically clearer. However, the nature of transition estimation 
has not changed; only the perspective and the tools have improved. 

At this stage of developments the concepts of chaos and strange 
attractors show little if any promise in this task. The name of our 
technological game is "convected spatial disorder in x- growing shear 
layers" not "temporal disorder in recycled closed systems of constant 
scale". Transition and turbulence in open-flow systems are 
evolutionary processes, not balanced, equilibrated states . The 
serious difficulties in applying these concepts to high-Re open flows 
should not detract from their great conceptual and philosophical 
significance in less complex, more cohesive nonlinear systems across 
scientific disciplines. 

It is impossible to trace the myriad of discussions, letters and 
papers which have influenced the perspective offered in this review. 

To all those contributors the author is grateful for keeping him 
curious and young in retirement. Many are acknowledged indirectly by 
citation of their papers. The author also wishes to acknowledge 
partial financial support in preparation of this review under AFOSR 
Grant No. 86-0165 and a grant from ICASE and the Computational Methods 
Branch of NASA Langley Research Center. 
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APPENDIX: A TASTE OF CHAOTIC BEHAVIOR 

(A-l) New concepts and philosophi cal implications 

Chaos theory uses language far removed from engineering 
terminology. The most accessible exposition of the field is found in 
the books of F.C. Moon (1987) and Berge - ", Pomeau and Vidal (1986). Here 
we propose merely to illustrate how a nonlinear Van der Pol oscillator 
can undergo chaotic motion in response to pure harmonic forcing, a 
behavior seldom discussed in engineering courses. The abstractions of 
chaos may not impact practical engineering concerns in turbulent 
systems for a long time. Conceptually, however, they provide a new 
basic insight: erratic or chaotic behavior is to be expected across 

disciplines in DETERMINISTIC nonlinear systems with just a few degrees 
of freedom even in absence of random disturbances . 

Turbulence was discovered more than a century ago and fluids 
engineers have known over forty years that its onset occurs after just 
a few instabilities. The recognition that qualitatively similar 
aperiodic behavior takes place commonly in sufficiently nonlinear 
mechanical (e.g. vibrations), electrical, biological, chemical, 
ecological, etc systems provides a unifying philosophical perspective 
on our world. That fibrillation in the human heart heralds the onset 
of a chaotic mode of this electromechanical system represented a 
revelation which in turn helped the design of defibrillators. Darwin 
too would have found the concepts of chaos stimulating. 

( A-2^ Limit cycle - attractor 

The Van der Pol oscillator mimics flow instabilities in allowing 
small displacements y at small velocities y draw energy from the 
system and grow by virtue of the locally negative damping coefficient 

2 * 

c ( y -1) of y in the equation shown in Fig A-l . In the unforced case 
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(1), b - o, the motion is thus self -excited . For this second-order 
system, we need to specify the initial values of y and y (termed phase 
or state variables) to evaluate the constants of integration and obtain 
a specific solution in time. As Fie A-lb shows, an initially small 
displacement amplifies rapidly for the parameter value c - 1 before it 
"saturates" at a fixed amplitude. In an autonomous system (one without 
explicit time dependence in the equations) it is possible to express 
one phase variable in terms of the others at all times and plot them 
against each other, as in Fig A-la for the 2 -phase space. 

The specific time -trace solution of Fie A-lb corresponds to the 
inner trajectory, unwinding from the origin in Fig A-la toward the 
saturation curve C, called the limit cycle . For any initial state of 
the oscillator, i.e. any point y(o), y(o) within the contour C, there 
is a unique inner trajectory through this point. Each 
solution- trajectory unwinds towards larger y and y values and 
approaches C asymptotically. An infinity of these trajectories 
corresponding to all initial state variables inside the contour, crowds 
outward toward C from within. Solutions associated with initial y(o) , 
y(o) values outside of C have positive damping, lose energy and wind 
down towards C, the common asymptote. Only two are shown in Fig A-la . 
but again there is an infinity of these outer trajectories crowding 
down towards C. Given the y dependence of the damping, the geometry 
and direction of the trajectories is set. In a 2 -phase space the 
long-time solution- trajectories are thus constrained to approach a 
single finite -amplitude closed curve, the periodic limit cycle . the 
dominant intrinsic characteristic of the nonlinear oscillator. These 
long-time solutions are insensitive to initial conditions. The 
Karman- &enard vortex street close behind a cylinder for 45 < Re < 160 
exemplifies such a periodic attractor in an open- flow system. 

(A-3) Oscillator forcing - threedimensional phase space 

When the Van de Pol oscillator is forced harmonically, the term 
rbc«>coswt on the right side of the equation in Fig A- 1 renders the 
system non- autonomous . To obtain solutions- trajecteries with 
properties similar to those in the 2 -phase space of Fig A-la . we 
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eliminate the explicit t-dependence of coswt by introducing a third 
phase space variable through 9 - u>. The forcing term becomes cos0 as 
can be seen in the equivalent set of three first-order equations in Fig 
A- 1 . Again, for any combination of the parameters £, b and cj, there is 
a unique traj ectory- solution through every point of the 3 -phase space 
determined by the prescribed initial values x, y, and 6 at t - o. 

For some parameter combinations, as the transients subside the 
solutions are asymptotically attracted to the toroidal surface in Fig 
A-lc - the generalization of the one -dimensional limit cycle. The 
quasiperiodic laminar wake oscillations in Fig 12b correspond to such 
an attractor with two noncommensurate frequencies f^ and see S tuber 
and Gharib (1988) for more details. 

For some parameter combinations, however, the trajectories become 
disorderly, as if drunk with the freedom of motion provided by the 
extra phase dimension. The contorted family of trajectories still 
fills the 3 space, but is locally more dense in some regions, as it was 
in the proximity of the limit curve C in Fig A-lb . As time increases 
indefinitely the transients die out and the trajectories approach an 
attractor subsnace not describable in terms of classical geometry. 

This " strange attractor " is made of an infinite discrete set off 
convoluted, nearly parallel surfaces with spacing becoming locally 
infinitesimal . 

The asymptotic set does not fill the 3 space completely; nor is it 
a finite sum of two-dimensional surfaces. The attractor has an 
effective fractional dimension between 2 and 3. Various determinations 
of the "fractal dimension" of a strange attractor naturally require 
limiting procedures and large amounts of numerical or empirical data, 
with attendent problems in accuracy, Moon (1987). Yet the attractor 
dimension is a quantitative intrinsic property of the given nonlinear 
system which distinguishes it from other systems. A fluid system 
governed by Navier- Stokes equations could potentially have very large 
dimensions so that some critics question inferences of behavior based 
on experiences with low- dimensional systems. Sreenivasan (1985) 
estimated dimensions on the order of 20 in a cylinder wake at Reynolds 
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numbers near 10,000. 

In forced chaotic systems, sampling of the state variables at 
discrete periods 2nn/<j) corresponds to so called Poincare' sections 
labeled P in Fig A-lc through which one can study the evolution of the 
system in more detail. In this figure the signature of the 
quasiperiodic regime i.e. intersections of the asymptotic trajectories 
with the Poincare cut, are seen to aggregate into the closed curve Q. 
For the chaotic regime, Figure 3-14 of Moon (1987) shows the 
intersections of the asymptotic trajectories to form an infinite set of 
"highly organized points arranged in what appear to be contorted 
parallel lines". When one enlarges a portion of the section, similar, 
smaller-scale structure is revealed, a characteristic of fractal sets. 
This phase-space structure is reflected in the aperiodicitv of time 
traces of any state variable and their broad-band spectra . An 
oscillator, driven sinusoidally at constant frequency may respond in 
gibberish, not noise -induced, but deterministic gibberish! The Duffing 
oscillator, with linear positive damping and a cubic spring element 
responds similarly. 

The equations are relatively simple; even the crucial nonlinearity 
in the Duffing case is simple. Yet predictability and reconstruction 
of the past is thwarted . In principle, the results are repeatable. 

But there exists extreme sensitivity to initial conditions which 
control the specific solution- trajectory. Neighboring trajectories in 
chaotic systems tend to diverge exponentially in time, the divergence 
being characterized by positive Lvapunov exponents - another 
"measurable" property of the system. A minute change in the initial 
value, say of the velocity y, may therefore cause a shift to a 
different region of the attractor and different asymptotic gibberish. 

(A-4) On implications for open- flow systems 

The nonlinear behavior of systems with two and three degrees of 
freedom, summarized in Fig A-l . departs significantly from that 
stressed in undergraduate engineering texts . On one hand there is the 
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powerful constraint of the limit cycle and toroidal attractors with 
their lock-in, threshold, jump and hysteresis phenomena, elements of 
which are evident in all fluid oscillators. On the other hand, there 
is the intrinsic chaotic behavior which may be a low- dimensional 
manifestation of phenomena somehow related to turbulence in the 
high- dimensional Navier-Stokes flows. We first note that Fig 12c 
demonstrates the existence of chaotic wake motion without 
threedimensional turbulence. If the motion in that laminar wake could 
be maintained twodimensional, it would be difficult to distinguish it 
from theoretical twodimensional turbulence (with possible inverse 
cascade processes) - if a distinction does exist. Such a 
twodimensional chaotic flow would be a realization of a Navier-Stokes 
solution satisfying the turbulent syndromes of Fig 3 except that 
vorticity would remain two dimensional . The solution would be 
determined by vorticity dynamics without the powerful tilt-and-stretch 
mechanism of threedimensional turbulence. The common turbulent wake of 
Fig 12d would presumably have a different, higher- dimensional 
attractor. But how to translate the attractor characterization into 
information in the physical space, of general or specific engineering 
interest such as entrainment and mixing in the two wakes above? By 
building up a dictionary of strange attractors? 

A major problem is that the theory of dynamical systems has 
concentrated on evolution in time, and by virtue of its tools has 
nothing to say about the structure of chaotic flow systems in the 
physical space. The paradigm is one of temporal chaos. Yet in 
shear- flow turbulence threedimensional spatiotemporal organization 
determines the varied transport properties which are of primary 
interest in engineering. The limited experimental confiramtion of 
chaos theories in flows, Guckenheimer (1986), comes from closed-flow 
systems: flows between rotating cylinders and the Benard- Rayleigh 

layers heated from below. In both, the vorticity is recirculated upon 
itself ("working upon itself"). Streamwise change of scale and 
streamwise convection of vortical structures do not play a role as in 
open- flow systems. 


- 57 - 


A-6 


In important open-flow systems vortical structures move downstream 
while gradually spreading across streamlines. Besides this parabolic 
behavior in the mean, there generally exists streamwise nonhomogeneity 
and associated intermittency in the structures comparable to the shear 
layer in scale. Vortical structures are not recycled; there are no 
associated iterative processes which are basic to the theoretical 
generation of strange attractors. In a Benard-Rayleigh cell or the 
fibrillating heart, the same object is interacting iteratively with 
itself. In open flows the objects are passing by. Consequently, if 
strange attractors can be defined for x-dependent systems, they are 
likely to be glimpsed by an observer traveling with the objects. 

R.J. Deissler (1985) concluded on the basis of numerical 
experiments with the nonlinear Ginzburg -Landau partial differential 
equations in x and t with a convective term, UdA/dx, that if Lyapunov 
exponents were to be meaningful, the underlying measurements should be 
carried out in the frame of reference moving with U. That measurements 
by fixed probes in a physical shear layer sense primarily the effect of 
nonhomogeneities from upstream rather than some intrinsic properties of 
the vortical entities has been recognized among engineers since Favre 
and coworkers first measured space -time correlations in turbulent 
boundary layers in 1953. 

An associated issue is that of external noise. Noise in flow 
systems is omnipresent, broad-band and hard to distinguish from 
manifestations of an intrinsic strange attractor (except when kept very 
low). However, noisy fields, nonhomo gene ous in x, can be selectively 
amplified by open- flow laminar shear layers and can give rise to strong 
nonlinear fluctuations sensed at any x, Sections 2 and 3. Furthermore, 
we have seen in Fi g A-l that external harmonic forcing of sufficient 
magnitude can actually modify the dimensions of the system and bring 
about chaos where there had been stable periodic behavior. In 
open- flow systems small -amplitude noise can thus lead to large randomly 
modulated fluctuations which can add to and even modify the character 
of the putative intrinsic strange attractor of the shear layer. 
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This is not the case in closed flow systems where noise generally 
adds a low-level broad-band background, except when several modes of 
incipient instability compete at a critical stage and perhaps in the 
more erratic cases of large aspect ratios, Guckenheimer (1986), 

Deissler and Kaneko (1985) mimicked the open-flow noise amplification 
with the cited Ginsburg- Landau model even without the troublesome 
x-dependence of the shear-layer thickness. They conclude that the 
"noise sustained structures" are both seeded by the noise as well as 
destabilized by the remaining noisy features after saturation. While 
the secondary instability in their model is rather specialized 
(side-band instability) their generalizations agree with our account of 
primary and secondary instabilities and their inductions in Sections 2 
and 3.3. 

The issues raised in this Section prepare the discussion in 
Sections 5.2 - 5.5 where they are further elaborated in context of 
specific families of flow fields. 
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INSTABILITIES -+- TURBULENCE 


Dynamics & Kinematcs of VORTICITY DISTRIBUTIONS 

(BIOT-SAVART IMAGERY) w(x,y,z,t)=j*rjn-kf 


U 

§ 




initially quasi-homogeneous in x and z 


SEQUENCES OF 

some simultaneous 


INSTABILITIES 
co RESTRUCTURING 


1,2,3... 


First 

Primary 


Instability 


in x 


MECHANISM: inflectional 


viscosity 


tuned 


TURBULENT 

STRUCTURES 


= Dehomogenization — Restructuring 

._X 

conditioned 


in z 
centrifugal 


- X= . 

°r 

w=kAf 


Yc, ^d 

iTIT 


-C _ c -• 

fini mi i 7u»iHiiiw»wn 


wskAf 



u. =i A( 


(FIRST) NET VORTICITY CHANGES AVARAGE TO ZERO ZAco-0 


Initially 

linearizable 


cosa r (x-c r t)e CrC|, 


-a,x 


cos« f (x-ct)e~ 1 


cosflze' 


a,x 


followed by slower nonlinear growth 

Ns ^possible approach to ‘saturation’ 

ss'neutral’ max. of first restructuring 


USUAL ASSUMPTION (not sufficiently general) 


f f 


secondary Instability from local nonlinear equilibrium 

initially linearizable around this 


X 


linear D.E. with periodic coefficients 


FLOQUET SYSTEMS 5 PARAMETRIC INSTABILITIES 


FIG. 2 VORTICITY RESTRUCTURING IN INSTABILITIES 
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Examples 


SECONDARY : pairing in mixing layers 


/ multiple amalgamation 
l when stimulated by f/n 


tilt & stretching in C.H&K breakdowns in BL’s 
rotating wavy Taylor cells in Couette fijflow 


TERTIARY : hairpin-vortex formation (3D inflectional ) in K breakdown 
I 2-frequency wavy Taylor cells 


I 


y ULTIMATELY OPERATIONALLY INDISTINGUISHABLE FROM 
TURBULENT STATES characterised by 4 key SYNDROMES 


© irregularity (disorder, partial) -^stochastic nonlinear systems 

3 Dim vorticity (eddying) -^-ensemble of loosely coupled 

deformable gyroscopes 

@ diffusion far in excess of molecuxar mixing -s».both friend and foe 


U on a large range of scales 
broad-band spectra 


(3) dissipation 


see-Tennekes &Lumley (1972) 

and film by R. W. Stewart (1968) and 8-page summary 

RESTRUCTURING CONTINUES in turbulent shear flows 

--s» large - scale coherent structures.fine-scale intermittency. 


© represents strange behavior in continuum mechanics - 

-apparently common in NONLINEAR DISSIPATIVE SYSTEMS with 
'strange attractors' larger number of degrees of freed 


om) 


EXISTENCE of Re m j nt for self sustaining wall turbulence 


on scales 


upstream of Re m j nt 

buffeted laminar BL's 
lacking true turbulent 
transport (3) 


corresponding to the thickness 6 of the layer 

presumably through intermittent *BURSTs" 
of the viscous sublayer, correlated disturbance 


motions producing 
negative Reynolds 
stresses-puv 


SWEEP 


EJECTION 
_ v+ (large) 

P u ' 

IIIIIII H IIIIIIII I II H IIIIII I llllllllllllllllllllllllllll 

mushrooms or 
llke hairpin loops 


IN HIGLY ACCELERATED TBL BL's 

the viscous sublayer is stabilized, bursts stop, 
the sublayer becomes a ‘buffeted lam. BL\ decoupled from the 
wake-like, decaying turbulent outer layer (STERNBERG 1954). 


This is one mode of RELAMINAR1ZATION , 

see NARASIMHA & SREENIVASAN (1979) for other modes. 


FIG. 3 HIGHER RESTRUCTURING AND TURBULENCE 
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FlOt plate Re ©mint ~ Re ©cr 200 

^ — <g) — 


Attachment line on swept wing 



Pipe 


i 


VT) 

J J 


Re ~ 2300 equilibrium puffs « Re 

Dmint __ 2 700 wall bursts, slugs C ‘ 


Large aspect ratio ducts 


2h 


Re nmint ~ 1500 non-decaying tbl patch «Re hcr - 5700 


Fig. 4 : Minimum Re for self-sustained turbulence compared to linear 

critical Re for T-S like instabilities. 
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LINEARIZED RESPONSE to 

LOCAL IMPULSE DISTURBANCE at x=0, t=0 


proportional to e 


j(kx-wt) 



CONVECTIVE INSTABILITY 

e.g. Boundary Layers 

GASTER SPOTS (1975) 


i 



ABSOLUTE INSTABILITY 

Onset Criterion: 
group velocity 
Cg=da>/dk = 0 
(double root) 


STABILITY u),<0 for all x,t: 

impulse response decays everywhere 


See HUERRE + MONKEWITZ (1985) JFM 1_5§, 1 5 1 

TEMPORAL — - linear response to continuous Asin«x at t=0 
llNSTABILITYl — dependent on INITIAL CONDITIONS in x and/or t 
SPATIAL — linear response to continuous Asin»t at x=0 


FIGURE 5 
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Boundary layer 

• convective' 


LOCAL LINEAR 
INSTABILITIES 

•ABSOLUTE 


V downstream 

PROPAGATION^ 

' some upstream: 

potential resonance in 

Bluff near wake 


x 


LOCAL means: 

at each x, replace x-varying shear layej by parallel, 
constant layer with local mean profile U(y) - neglects strong 
destabilizing adverse pressure gradient near body! 


for Bluff wake 

y 



Near: absolutely unstable 


y 

— , 

U(y) 


| 

past closure of 
recirculation pocket 





Farther: convectively unstable 


RIGOR in usage of local mean profile of vigorously oscillating wake? 


PROBLEM of NONLINEAR UPSTREAM EFFECTS through PRESSURE : 


y 



This unsteady feedback converts the linearly convectively unstable 
mixing layer into a weakly (non-robust) absolutely unstable one! 


[GLOBAL INSTABILITY 


arises from 'resonance’ in x when 


non-parallel x evolution is directly taken into account. 


Modeling via Ginsburg-Landau equation, CHOMAZ, HUERRE 
REDEKOPP (1988, Phys. Rev. Letters, 60, 25) show: 

NEED long enough x range of Absolute Local Instability to build up an 
unstable eigenfunction in x on top of the local y eigenfunction of 
the layer. 

PROBLEM of Boundary Conditions in x, especially at body 


GROWTH TO NONLINEAR SATURATION and ENERGY BALANCE 


SYSTEM then ROBUST to environmental disturbances! 

HYDROACOUSTIC GLOBALLY UNSTABLE SYSTEMS 


e.g. jet edge-tone: 

2^2**. 

Dt — » U 




pressure 

,<lipole 

~ ~ ~ ~ ^'feedback 
L * to lips 



IRREGULAR f JUMPS 


L 


FIG. 6 CONVECTIVE, ABSOLUTE and GLOBAL INSTABILITIES 




MIXING LAYER INTERACTIONS 
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* * 



FIGURE 7 Summary of primary and secondary instabilities and their upstream influence on artx- J 
ficially unstimulated disturbances at origin of the layer; conditions correspond to the inset, 
courtesy of Drubka (1982) . Locking of subharmonic f/2 onto the nonlinear fundamental f at A* 

initiates the secondary instability. TS designates unstable modes of the primary instability. 
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INSTABILITIES AND ROAD TO TURBULENCE in 

OPEN-FLOW SYSTEMS 

ARE FORCED by ENVIRONMENTAL DISTURBANCES 

Cexept in rare cases of global-absolute instabilities) 

R E_CE_P T I V I T Y ; conversion of external disturbances into 

internal homogeneous growing vorticity modes 


SMALL DISTURBANCES 


LARGE DISTURBANCES 


Linearizable Receptivities: 

• in principle, SOLUTIONS to 


NONHOMOGENEOUS 


H Orr-Som. / 
Gortler > 
cross-flow) 


SYSTEM: 


PARTICULAR (FORCED) SOLUTION 

PLUS 


all HOMOGENEOUS SOLUTIONS 
needed to satisfy together all 
BOUNDARY CONDITIONS in space 
and time, EXACTLY (not on the 
avarage, in case of turbulent 
disturbances) These conditions fix 
the coeff C« in the symbolic form: 

S °W + JCu.HomoSol,, 

For larger x and t all but one 
(or few:j) HomoSolp and SoWt 
decay, leaving the dominant 
response Gj .HomoSolj 

Ex: BECHERT (1988) JFM 186.47 


Concept of NONLINEAR Re m j„t 
for self-sustenance of turbulence 
in a given Boundary Layer, 
generally below linear Re« 

UNKNOWN NEW 

x.y.z.t - dependent 

BASE FLOWS 

If locally linearizable: 

RAPID INFLECTIONAL INST ABILITIE! 
(with forcing) 


BYPASSES 


of previously explored linearized 
fields and mechanisms. 

EXAMPLES 

• Large, single-PROTUBERANCE 
instability and transition 

• Impingement of wakes from 
upstream in turbomachinery 

• Turbulent patches in 2D 
Poiseuille ducts 

• Larger distributed ROUGHNESS 


see Morko vin( 1 984) ‘Bypass 
Transition’, NASA CP-2386, 
Sympo: Transition in Turbines 
NASA Lewis Res. center 




3ASS: Enhanced TS and Herbert receptivities to tree-stream fields 
due to mean-flow changes BELOW the heights of smaller distributed roughness, 
CORKE, BAR-SEVER , MORKOVIN( 1 986) Phys. Fluids 29,3199 


FIG. 11 ENVIRONMENTAL FORCING & RECEPTIVITY 
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11.5 cm/sec 
in water 





unsteady vorticity 
forced at frequency 
by heating elements 


(V • : PAGE IS 

OF POC ALITY 


fi locked on 
natural-street 
frequency f„ 

= 4.075 Hz 




chaotic laminar 
f 1 =3.55 Hz 
f 2 = 4.075= f„ 
f 3 =5.42 Hz 
some 3D motion 



no forcing ; 
turbulence tripped 
symmetrically 
by wires in B.L. 

U = 1 3 cm/sec 



COURTESY: K. Stuber and M. Gharib 


FIG. 12 EXCITED STATES IN SLENDER WAKE 
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SOME ROADS to TURBULENCE ± CHAOS 


Q: can shear layers growing in x have STEADY strange attractors ? 
Q: can chaos theory accomodate lam.-Tu coexistence? 


U 


45< Re<160 



Re>200 
Tu downstream 

* 

marches upstream 


GLOBALLY unstable system NEAR BODY 
' ' some local feedback through unsteady backflow, 
stable, LOCKED ’clockwork’ oscillator 
< 3D-instab. (absolute?); role of wi disturbances? 

| sym.-antisym. mode competition? GerrardC 1 978) 



Tu reaches fluid in first shedding vortex 

Tu in separating layer causes c 

* reattachment and bubble; Re~2.l0° 

^ IMPORTANT SPANWISE NONUNIFORMITY 

Tu ^onset upstream of lam.sep., y 

crooked tbl.sep. line, low C' L ; Re~1.5.lO D 

Tu marches upstream in thin B.L.; 

roughness sensitive; wake widens, 
mild C' L rise as sep. line straightens 


bubble 




B.L.: CONVECTIVELY unstable RESPONSE to low free-stream Tu: 

superposition of lam. TS spots, amplifying quasi2D in x 



As U increases, pattern moves upstream, 
but Re tf not constant! Unit-Re effects ! 


UNLIKE for pure STRANGE ATTRACTORS. 

SIGNAL IRREGULARITY strongly NOISE DEPENDENT 

Q: Any way of experimentally sensing the intrinsic SA irregularity if present? 
Q: Invariance of such intrinsic measures in y and z, and especially in x? 


FIG. 13 EVOLUTION TO TURBULENCE IN 

GLOBALLY AND CONVECTIVELY UNSTABLE EXAMPLES 
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of wave packets.- 
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components of 
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IRREGULAR TERTIARY 
INSTABILITIES 
e.g 3D inflectional- 
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VORTICES 
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uncharted interactions 
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FIGURE 15. 1984 system portrait of environmentally stimulated 
vorticity perturbations on their evolutionary paths to 
turbulence in undistorted boundary layers and ducts. 
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mm 


Forced Van der Pol oscillator: 

y + « (y^ - 1 )y + y = b<ucosui 


C 1 J unforced case b=*0 

autonomous system 


OUTERi 

rTRA JECTORIES in phase space 
inner 'both attracted to the same 

set of points: LIMIT CYCLE C 


therefore motion becomes 
periodic in time , has 

TIME RECORD 


(2) forced case b*0 

To make system autonomous (free of explicite t-dependence) 
introduce a third dimension (phase variable) 0 by setting G = a>. 
Equivalent system of three first-order equations is: 

y = e | x - ( y®/ 3 - y ) + bsin0| ; x =-y/« ; 0= u> : PARAMETERS «,b,w 


Solutions trajectories attracted to: 

curves on a limiting TOROIDAL SURFACE—* 
STRANGE ATTRACTORS: collection of 
infinite sets of convoluted ‘parallel" surfaces 
with spacing locally approaching infinitesimal, 
studied by Poincare’ sections, P—» 

PROPERTIES have aperiodic , bounded traces 
in time, corresponding to broad-band spectra 
Positive Lyapunov exponents trajectory 
divergence 

Supersensitivity to initial conditions in phase 
space 


After Me Lachlan (1947), Oxford Un. Press and Berge’ 
Pomeau and Vidal (1986), John Wiley & Sons 



FIG. A- 1 ROAD TO CHAOS IN A FORCED 

NONLINEAR OSCILLATOR 
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